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The history of mankind is one of continuous development from 
the realm of necessity to the realm of freedom. This process 
is never-ending. In any society in which classes exist class 
struggle will never end. In classless society the struggle 
between the new and the old and between truth and falsehood 
will never end. In the fields of the struggle for production 
and scientific experiments mankind makes constant progress and 
nature undergoes constant change j they never remain at the 
same level. Therefore , man has constantly to sum up experience 
and go on discovering , inventing, creating and advancing . Ideas 
of stagnation , pessimism, inertia and complacency are all wrong. 
They are wrong because they agree neither with the historical 
facts of social development over the past million years, nor 
with the historical facts of nature so far known to us ( i.e . 
nature as revealed in the history of celestial bodies, the 
earth, life, and other natural phenomena ). 
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STROPS IS 


The fission process continues to be a fascinating and 
challenging area cf study even though almost four decades 
have elapsed since its discovery. During this period a 
large number of investigations on both the theoretical and 
experimental aspects of this process have been carried out. 
However, a clear picture, coherently linking all aspects of 
this complex process, is yet to emerge. 

The charged liquid drop model corrected for shell 
effects has been quite successfully used in describing the 
main features of the fissioning nucleus, eg. the fission 
barrier with respect to deformation, fission isomers, 
intermediate structure in subthreshold fission cross sections 
etc. One of the crucial questions in the passage of the 
nucleus from saddle point to scission is the extent to which 
there is coupling between the collective and particle degrees 
of freedom. In other words, if the process is adiabatic with 
respect to the particle degrees of freedom, then the decrease 
in potential energy from saddle to s cas sxono -remains associated 
with the collective degrees of freedom and appears primarily 
as kinetic energy of the nascent fragments. On the other 
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hand, if the motion is non-adiahatic then there will be 
transfer of collective energy into nucleonic excitation 
energy, A reasonable explanation based on this kind of 
motion is provided by the statistical theory of fission. 

Present experimental evidence does not allow an unequivocal 
decision in favour of either the adiabatic or non-adiabatic 
model. More inf oimation on the properties of the nucleus 
at the saddle point and scission point is necessary in 
order to shed more light on this question. 

Experimental investigations on light-charged-particle- 
accompained (LCP) fission, also referred to as ternary fission, 
have indicated that, because of the time and place of 1CP 
emission; a study of this mode of fission would contribute 
significantly to the understanding of the scission stage 
of fission. 4- 

The present thesis consists of two parts. The first 
describes some experimental work examining the effect of the 
saddle point configuration on ternary alpha (IRA) yield and 
average energy. The second part describes computational 
work aimed at studying the scission configuration of the 
fissioning nucleus. 

The introductory chapter describes the salient 
features of the fission process. The characteristics of 
binary and ternary fission are compared in order to show the 
essential similarity of these two modes of fission. The 
emphasis on study of ternary fission is, however, motivated 
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by the fact that the LCP serves as a powerful probe of the 
scission process . j 

^Chapter II summarizes the channel theory of fission, 
developed by A. Bohr, and its application to the explanation of 
some important characteristics of low energy neutron induced 
fission. The high light of this model is its emphasis on the 
role played by the saddle states (or channels) in determining 
the post-scission characteristics of fission. The question 
of validity of the adiabatic model is intimately linked to the 
success of this theory in the explanation of such fission 
characteristics . 

The details of the experimental setup used in the LCP 

yield and energy measurements in low energy (thermal to 1 MeV) 

235 2^9 

neutron induced fission of U and Pu are described in 
Chapter III. Details of data analysis, specifically the 
corrections for chance coincidences and energy loss effects, 

«* 

are also given. 

The results of these measurements and their possible 

interpretation on the basis of the channel theory are 

presented in Chapter IV. The results on IRA yield and 

average energies and widths as a function of neutron energy 

indicate definite structure outside the limits of error. The 

235 

structure in yield in U is, however, more pronounced than 
239 

in Pu. The yield of the low energy component, attributed 
to tritons, also shows significant structure in the neutron 
energy region investigated. The results on LEA yield are 



xviii 


discussed in terms of the spins, /parities and fission cross 
sections of the low lying saddle point states of these 
nuclei. However, no definite conclusions regarding the 
'triton' yield could he "drawn because of large errors due 
to poor statistics. 

The method of trajectory calculations has been used 

to study the scission configuration. The post-scission 

motion of the fragments and alpha particle were simulated on 

a computer and by an interpolation method the scission 

configuration was obtained. Chapter V describes the method 

\ised and the calculations carried out for spontaneous 

252 

ternary fission of Cf. This nucleus was chosen because 
of the vast amount of experimental results available. The 
method used in this work differs from almost all previous 
methods in that no apriori restrictions were Imposed on the 
initial parameters used as the starting point of the simula- 
tion i.e. in the approximate description of the scission 
configuration. The results indicate that such a method is 
valid and it has been possible to obtain hitherto unobtained 
distributions of initial parameters describing the scission 
configuration and some correlations between them. 

Chapter VI contains a brief summary of the results of 
the experimental work and their interpretation and the results 
on the scission configuration obtained from trajectory 
calculations. 



CHAPTER I 


INTRODUCTION 


1 . 1 The fission process 

The fission process, that is, the breaking-up of 

a nucleus into two or more heavy parts, was discovered about 

40 years ago by Hahn and Strassman (1938). It has since 

been extensively investigated both theoretically and 

experimentally. The first explanation for this process was 

given by Meitner and Frisch (1939) followed by a more 

quantitative description by Bohr and Wheeler (1939). In 

general nuclear reactions the question of whether the process 

is through a direct interaction or via compound nucleus 

formation is an important issue in itself. In contrast in 

this, fission is a clear case of a process via compound 

nucleus fomnation. This is because in fission we deal with 

collective nuclear motion and the characteristic time for 

—21 

such motion (vibration) , 5x10“ sec., is long compared to 

-21 

the time for a nucleon to cross the nucleus 0.3x10 sec. 
The compound nucleus in question may be formed by absorption 
of a particle or electromagnetic radiation - as in induced 
fissions or it may even fission from its ground state 
(spontaneous fission. Fig. 1.1). The compound nucleus 
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formed by absorption of a particle would be in an excited 
state and the excitation energy is distributed among a large 
number of degrees of .freedom. The complex state of motion 
initiated thereby may be described in terms of collective 
nuclear rotations, vibrations and single particle excita- 
tions. In the compound system there is competition between 
the various modes of dissipation of the available energy^ 
neutron emission, radiation or division. For fission to 
take place, it is necessary that a sufficient amount of 
energy becomes concentrated in the potential energy of 
deformation to enable the nucleus to pass through the saddle 
point shape - a shape corresponding to maximum potential 
energy of deformation. At the saddle point the repulsive 
CSoulomb forces balance the attractive surface tension force 
(arising out of the nuclear interaction). Beyond the saddle 
point and till the scission point the repulsive forces 
increasingly dominate the attractive ones. The scission 
point is the point at which the nucleus divides into two or 
more smaller nuclei such that their nuclear forces cease to 
play a part as far as their mutual motions are concerned. 

Thus the fission process may be described as the 

succession of three different phases. The compound nucleus 

undergoes a long series of oscillations until one of them 

leads to the passing of the saddle point, followed soon 

—21 

afterward (<*'«' 10"” sec.) by the splitting-up of the nucleus. 
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In most cases this division is mass asymmetric leading to 

light and heavy fragments. These primary deformed fragments 

move in opposite directions under the influence of their mutual 

Coulomb repulsion which results in their kinetic energy. They 

de-excite by emitting 'prompt' neutrons and gamma rays in time 

—11 

intervals of 10 and 10 sec., after scission, respectively. 
[An analysis of the angular distribution of neutrons emitted 
in fission (Milton and Fraser, 1965) indicated that a small 
fraction of neutrons (which are isotropic in the laboratory 
system) are probably emitted at scission]. These fission 
fragments are neutron rich and far from the stability curve 
and decay by delayed neutron and gamma emission and by beta 
emission to stable end products. Most physical measurements 
are carried out on the fission fragments whereas radiochemical 
measurements are performed on the fission products. 

The total energy released in fission, comprising the 
kinetic energy of the fragments, the energy associated with 
neutron and gamma emission etc., is about 200 MeV and is equal 
to the difference in masses of the original compound nucleus 
and the' two fragments. The explanation of this energy 
release was first provided by Meitner and Frisch (1939) and 
confirmed experimentally by Frisch (1939). 

The most frequently occuring mode of fission is 
binary fission; in which the nucleus breaks up into two 
fragments. About 0.2^ of such events result in three 
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* Mission 

fragments - this division is called ternary fission, 
into four or more fragments is extremely rare. 

A number of experiments have been done to study 

■ rfrtw ott'pt' no theory , 

characteristics of the fission process. Ho^e > 
which comprehensively explains all aspects of this proce 
has emerged. Bohr and Wheeler (1939) developed a theory o 
fission based on the conception of the nucleus as a lil 
drop, Frenkel (1939) independently proposed a similar 
Their application of this theory accounted for a numbe 
features of the process, however, it did not explain the 
most striking feature of fission, namely, the asymmetric 

nature of the mass division. The liquid drop model ha 

. . . o-Hii usocL "today 

been developed (Swiatecki and others) and is 

as the starting point of calculations of binding energ' 

potential energy barriers (with respect to deformation) 

heavy nuclei. There are, however, marked deviations from 

the predictions of the liquid drop model on binding 

energies etc., in the vicinity of the well-known neutro 

proton shells. The first successful method for quantita 

vely introducing shell corrections to the liquid drop 

potential was introduced by Strutinsky (1967, 1968). 

procedure is now widely used in calculating fission barri 

etc. For actinides, calculations of this type give fi 

barriers with two humps separated by a second minimum 

(Fig. l.l) which provides the possibility for very 



(1 (Deformation) 

.atic illustration::- of single-humped. (— ) 
,oubl e-hum pod (— — } fission bar ri era . Intr: 
iations in the first and second wells are ■ 
tinted' class I and class IX' states f respect 
nsic c ha nit el 3 at the two harriers are 'also' 
'trpted , , $he transition in the shape of.'th 
us 'as a function .of . deformation Is sohemat: 
sente'd' in. the u per' part .'.of the. figure '[A: 
nbosch. 'and Huiaenna # ' 19T5 ] . : .. i * . , ■ 
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states (called class II states) to exist in the second well. 
The existence of such states can explain very well some 
important experimental results such as fission isomerism and 
structure in subthreshold fission cross-sections. 

1 . 2 Bi nary fission 

As mentioned earlier, the most commonly occuring 
division in fission is into two fragments of unequal mass. 

This asymmetric mass division is one of the most important 
characteristics of the process and is exhibited by all 
heavy nuclei with A> 229. The mass distribution is 
consequently double humped. Some nuclei, however, exhibit 
an anamalous mass distribution; Bismuth (in fact nuclei with 
Z483) exhibit predominantly symmetric mass distributions 
while Radium has a mass distribution which is triple humped 
indicating the presence of both symmetric and asymmetric 
divisions. A number of theoretical models have been suggested 
to explain these data, yet no one theory can explain all 
these observations. Other features of the mass distribution 
are that the mass distribution for a particular nucleus becomes 
more symmetric with increase in excitation energy and that the 
mass distribution depends on the mass number of the nucleus 
undergoing fission. The foimer observation has been explained 
in terms of the 'washing out' of shell effects at high 
excitation energies implying that the fissioning nucleus 
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’sees' only a liquid drop type of barrier (the liquid drop 
model predicts symmetric fission as the most probable mode). 

Another important characteristic is the division .of 
charge between the two fragments. Nuclear chemists have 
made many measurements to establish the character of this 
division of charge by measuring the yields of stable nuclei 
resulting from beta decay of the fission fragments. However, 
the quantity of most interest is the charge-to-mass ratio 
of the fragments just produced at scission prior to the 
emission of prompt neutrons i.e., of the primary fission 
fragments. The most probable charge (Zp) of the fragments 
have been determined by measuring the characteristic 
K-x-rays in coincidence with the fragments (Reisdorf et al. , 
1971 ). The K-x-rays arise from internal conversion 
associated with the deexcitation of the primary fission 
fragment s(lO~ sec., after scission). It has been observed 

that the -charge division is such that the lighter fragment 
receives a larger fraction of the available charge than 
would be expected if the charge-to-mass ratio of the 
fissioning nucleus were preserved in the fragments. Various 
postulates have been used to explain these observations eg. 
the equal charge displacement postulate (Glendenin et al. , 
1951 ) which can be expressed as (Zp-Z A ) ligli ^. = heavy 

where Z^ represents the charge of the most stable isobar 
of that mass chain. This implies that there would be fewer 
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beta decays to stability in the case of the light fragment 

(relatively proton rich) than in the heavy fragment. 

Another postulate is that the charge division is such that the 

potential energy at scission is a minimum and another which 

requires a maximum of the total fragment excitation energy 

at scission (Wahl et al. , 1962). One additional point to 

note its that the ratio of most probable charges of the 

252 

fragments is very close to the mass ratio, eg., for Cf the 
most probable mass and charge ratios are ^1. 31. 

The largest fraction of the energy released goes 

into the kinetic energy of the fragments, resulting from 

Coulomb repulsion. The fragment energy distributions are 

related approximately to their mass distribution through the 

relation M^/Mg = It is also noticed that the total 

fragment energy is a slowly increasing function of mass 

2 

number or more precisely of ‘ the parameter, Z / A. One 
of the characteristics of the kinetic energy distribution (in 
low energy fission) is the decrease for symmetric fission and 
a prominent peak for asymmetric fission leading to heavy 
fragments with A 5 * 132. This kinetic energy decrease has 
been qualitatively understood in terms of the relative 
increase in the number of prompt neutrons emitted in 
symmetric fission (this will be discussed later). The effect 
of increasing the compound nucleus excitation energy is to 
decrease the structure in this variation of kinetic energy 
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release. For nuclei lighter than Radium, which exhibit 
symmetric mass yield curves , the maximum kinetic energy 
release occurs for symmetric mass divisions and falls off 
with increasing mass asymmetry at a rate dependent on the 
compound nucleus excitation energy (Unik et al., 1969). 
However, in heavier nuclei, which exhibit asymmetric mass 
division, the total kinetic energy for the symmetric mode 
is lower than for the asymmetric mode. The dependence of 
total kinetic energy on compound nucleus excitation energy 
has been studied in detail and such studies indicate that 
the kinetic energy release in spontaneous fission is lower 
than that in thermal fission (Okolovich and Smirenkin, 1963) 
and for still higher excitations the variation in kinetic 
energy release is very small. This fact suggests • that a 
large fraction of the compound nucleus excitation energy 
goes into fragment excitation and not into kinetic energy 
implying that the extent of stretching of the neck (the 
Coulomb potential energy at scission) remains more or less 
the same as the compound nucleus excitation energy increases. 
For symmetric mass division the drop in kinetic energy 
could, therefore, be related to the increase in fragment 
excitation leading to increased neutron emission. 

Most of the energy release not appearing as kinetic 
energy of the fragments is dissipated by prompt neutron 
emission. The number of neutrons emitted is a direct 
consequence of the amount of energy stored in the fragments 
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either as internal (excitation) energy or as deformation 
energy converted to excitation energy as .the fragments 
collapse to their equilibrium shapes. The neutron yield 
increases with increase in compound nucleus excitation 
energy which can he understood in terms of the weak 
dependence of kinetic energy on this excitation energy. As ' 
mentioned earlier the total number of neutrons increases 
for symmetric fission (Silton- andi Fr&se&v-I196l'$ Apalin et al., 
1964 ) with a corresponding decrease in kinetic energy. One 
of the most interesting aspects of neutronemission is the 
striking (’saw tooth') variation in neutron yield as a 
function of fragment mass. (Fraser and Milton, (1954) 

Apalin et al. , 1965; Terrell, 1965). This behaviour led 
Terrell (1962, 1965) to the conclusion that the neutron 
yields are closely related to the def ormabilities of the 
nascent fragments. Closed shell and near-closed-shell 
nuclei prefer spherical shapes and are resistant to defor- 
mation resulting in lower neutron yields. The 'saw-toothed' 
variation in neutron yield is fairly rapidly washed out with 
increasing excitation energy (Bishop et al. , 1970). 

The angular distribution of fission fragments has 
been studied in the excitation energy region up to several 
tens of MeV for many targets, for a variety of projectiles. 
(Gordon et al. , I960; Vandenbosch et al. , 1961; Bate et al., 
1965; Simmons and Henkel, I960). The most striking 
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features of the experimental observations are (a) the fission 
fragments have largest differential cross sections in the 
directions forward and backward along the beam (b) the 
anisotropies are largest for the heaviest projectiles and 
smallest for neutron and proton bombardment (c) the aniso- 
tropies are approximately the same for odd A targets as for 
even-even ‘targets and decrease with increasing Z /A. The 
angular distribution of the- fragments depends on the angular 
momentum brought in by the projectile and the fraction of 
this converted into orbital angular momentum between the 
fragments (characterized by K in Pig. 2.3). The anisotropy 
has been successfully explained using the channel theory 
(discussed latter) using the assumptions that the fragments 
separate along the symmetry axis and that K is a good quantum 
number beyond the saddle point. Such anisotropies have also 
been observed in photofission (Baerg et al. , 1959) however, 
the differential cross sections are largest in the directions 
perpendicular to the incident photon beam (electric dipole 
absorption) . 

Prompt y rays emitted from the fission fragments 
provide useful information about properties of the fragments 
at scission. There are typically about 8 photons per fission 
(in spontaneous fission it is about 10 ) with an average of 
1 MeV per photon. Both the average number (the multiplicity) 
and average energy per photon are functions of fragment mass. 
These y rays differ from y rays associated with neutron 



12 


capture because of the large amounts of angular momentum 
they carry away from the fragments ('-'10-12 h). Extensive 
investigations of the characteristics of these fission 
y-rays have been carried out (see for example the review in 
Vandenbosch and Huizenga , 1973) and have provided useful 
information about the de-excitation mechanism of the 
fragments. 

1 . 3 Te rnary fission 

The term 'ternary fission' will be used with regard 
to that division resulting in two heavy fragments and a light 
charged particle (1CP) first observed by Alvarez (Earwell 
et al. , 1947). Such a division normally occurs once in 
about every 500 fissions. Another possible mode of division 
is that resulting in three fragments of almost equal mass , 
however, such a division is very rare (Muga and Rice, 1969). 

A very large percentage ("'90%) of these LCPs are long range 
alpha particles (LRA) having a mean energy of about 15 MeV 
(Gosper et al. , 1967; Nadkarni and Kapoor, 1970), the 
remainder are mostly tritons, with an average energy of 
~8 MeV, and other particles such as protons, deutrons etc. 

The present study concentrates on certain aspects of. IRA 
fission. 

One of the most striking features about LCP 
emission is the angular distribution with respect to the 
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fragments. The fact that it is peaked in a direction 
perpendicular to the direction of the fragments (Titterton, 
1951) has led to the belief that the LCP is emitted very 
close to the moment of scission from the 'neck' connecting 
the two fragments. Once these LCPs are emitted their 
motion would be influenced by the Coulomb field of the 
fragments and their final energies and angular distribution 
would be sensitive to the separation and speed of the 
fragments at the instant of their emission. Therefore, a 
study of LCP emission is expected to provide useful infor- 
mation about the scission-stage of fission. 

The excitation energy dependence of LCP yield is 

fairly weak, decreasing by about 25 % for thermal neutron 

fission (E*~ 7 MeV) compared to spontaneous fission (E^O) 

e.g. , for 240 Pu (Nobles, 1962). At excitation energies over 

15 MeV there is some evidence of an increase in yield, 

(Thomas and Whetstone, 1966). It has also been noticed that 

the yield increases with increase in fissility parameter, 

2 

Z /A (Nobles, 1962). This is consistent with the idea that 
the yield is dependent on the amount of deformation at 
scission i.e., a ’neck' should exist; which, as indicated 

O 

by the liquid drop model, increases with increasing Z /A. 

The kinetic energies of the various LCP vary 
quite significantly (Table 1.1). In 235 U+n th fission the 
most probable LRA energy is 15 MeV (Nadkarni and Kapoor, 197C) 
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while the average triton energy is about 8 MeV. It has been 
observed that the average kinetic energy of the alpha 
partiole increases as the angle of emission changes from the 
most probable angle (see Pig. 5 . 5 ) • 

Another striking feature of LRA fission is the 

observed anti correlation between the average total fragment 

kinetic energy and the alpha energy. The value of this 

anticorrelation, averaged over all mass ratios, has been 

measured and found to be -0.44 i.e. = -0.44 

(Praenkel, 1967; G-azit et al. , 1970; Mehta et al. , 1975). 

252 

This anticorrelation decrease with mass ratio; in Cf spon- 
taneous fission it is about -0.2 in the far asymmetric 
region and increases to about -0.8 in the symmetric region 
(Mehta et al. , 1975). This experimental fact will be used 
later in connection with trajectory calculations used to 
study the scission configuration. 

1.4 Comparison of Binary and Ternary fi ssion 

A comparison of the mass -distribution in these two 

235 

modes of fission e.g., in thermal neutron fission of U 
(Schmitt et al. , 1962) shows that the low mass sides of both 
fragment peaks are identical within the error of measurement 
whereas the high mass sides of the peaks are displaced 
towards the low mass side in ternary fission. The net 
shift in the peaks amounts to about 4 mass units corresponding 
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to the mass of the alpha particle. 

The average total kinetic energy of the fragments 
in binary fission is almost equal to the sum of the average 
kinetic energies of the fragments and alpha particle. The 
average number of neutrons emitted in IRA fission of 2 ^ 2 Cf 
is 0.65 neutron less than in binary fission (STardi et aJ, 1970 
Piekartz et al. , 1970). Using this figure it has been 
estimated that the excitation energy of the fragments is 
5.7 MeV less in LEA fission than in binary fission. A 
remarkable similarity in these two modes is the 'saw tooth' 
variation of average neutron yield as a function of 
fragment mass. The total y-ray energy in binary and ternary 
fission appears to be the same to within less than 1 MeV 
(Adamov et al. , 1967? Ajitanand, 1969). These facts go to 
show the essential similarity of these two modes of fission. 

1.5 T he LEA as a probe of the fission process 

There have been various conjectures about the 
extent to which the properties of the nucleus at the saddle 
point influence the various fission characteristics. Por 
eg. the occurrence of negative parity states at the saddle 
point, for symmetric deformations, has been related to 
asymmetric fission. Another theory - the statistical 
theory, concentrates on the nucleus just at the moment of 
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scission where the characteristics of the process are 
believed to be decided. There does not exist adequate 
experimental evidence to decide to what extent the saddle 
states or statistical factors (at scission) influence the 
fission characteristics. 

The present thesis reports experimental results on 

the yield and average energy of LRA in low energy neutron 

induced fission of U and Pu. The variation of IRA 

yield is sought to be linked up with some properties of the 

nucleus at the saddle point, thereby providing some more 

information which could help in moving toward a decision on 

the above questions. By using the method of trajectory 

calculations, the scission configuration in ternary fission 
252 

of Cf has been investigated. The results on the parameters 

describing the scission configuration, such as energies of 

the fragments and alpha particle, interfragment separation 
* 

etc. , could also throw some light on such questions. 



CHAPTER II 


CHANNEL EFFECTS I N F ISSION 

2.1. In troduction 

A complete treatment of the fission process requires 
a dynamical study of the penetration of a multidimensional 
fission barrier. However, theoretical calculations mostly 
concentrate on the penetration of the potential energy 
barrier with respect to a deformation parameter. The poten- 
tial energy surface (or fission barrier) with respect to the 
deformation parameter cannot be obtained, at present, using 
microscopic calculations which take into account the detailed 
structure of the nucleus. This is because of a lack of 
precise knowledge of the effective nucleon-nucleon interac- 
tion in the nucleus and also because of computational diffi- 
culties both in complexity and length. In addition, one 
would be faced with the problem of defining the deformation 
coordinates. Nevertheless, calculations of the fission 
barrier using Hartree-Fock methods have been carried out 
(Brack and Quentin, 1973 j Flo card et al. , 1973); but despite 
their considerable interest and the rapid improvements in 
accuracy such methods cannot as yet provide realistic 
barrier shapes. 
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Certain macroscopic methods based on classical 
concepts have been used with remarkable success. The first 
such calculation was made by Bohr and Wheeler (1939) employ- 
ing the liquid -dr op-model (LEM) to study the energy systema- 
tics of the fission process. In this model the nucleus is 
represented by a charged incompressible liquid drop and its 
potential energy is calculated in terms of the energies 
associated with the volume, surface area. Coulomb repulsion 
etc. The latter terms are shape dependent. These terms are 
summarized in the Weizacker (1935) semi-empirical formula 
for the nuclear ground state energy t 

B = -D rj + O' + u A 2 / 3 + U (A-.?. 2 ;). 2 + 5(2 A) 

v o a 1/3 s T 4A 

where the first term is the volume energy tern, the second 
the Coulomb repulsion energy term, the third the surface 
energy term, the fourth a term depending on the neutron- 
proton symmetry factor and an empirically observed term 
depending on the odd or even character of the nucleon numbers. 
Prom this formula it is possible to show that the energy of 
a heavy nucleus is considerably greater ( /N '200 MeV) than 
the sum of the energies of two medium-weight nuclei of the 
same total mass, explaining the energy release in fission. 

It is also possible to obtain the ratio of the strengths 
of the repulsive coulomb forces to the strength of the 
restoring surface-tension forces for which the drop is in 
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unstable equilibrium - any small deformation will cause it 
to undergo fission. This ratio, denoted by the fissility 

p 

parameter Z /A, has a critical value of about 50. 

Much effort has been expended in studying the poten- 
tial energy surface of the deformed liquid drop. Most 
studies (Bohr and Wheeler, 1939; Frankel and Metropolis, 1947; 
Cohen and S’wiatecki, 1962) have employed the coefficients of 
a spherical harmonic expansion of the radius of the drop as a 
set of parameters describing deformation. If the potential 
energy surface is plotted as a function of one such deforma- 
tion parameter the curve shows a maximum .(Big, 1.1) and 

if it is plotted as contours on a surface using as coordinates 
the two principal parameters of symmetric deformation, Pg an< ^ 
p^ ; a saddle point would result for some ( p ^ ^ P4. ) (Lynn, 1968). 
A similar saddle point would be apparent in a contour plot on 
a many dimensional surface defined by coordinates which 
included also the higher order' coefficients of the expansion. 
It is the principal object of liquid -drop calculations to 
find the energy at the saddle point with respect to the 
energy of the undeformed drop, for this is the threshold of 
the fission reaction. 

Despite the qualitative success of the liquid-drop 
model in explaining the fission phenomena it could not 
account for may observed characteristics such as mass 
asymmetry, isomer fission etc. The correction of the 
liquid -drop model for shell effects provided a great impetus 
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to the study of fission. Such corrections* "were first 
successfully incorporated hy Strutinsky (1967, 1968) by an 
ingenious combination of the macroscopic and microscopic 
properties of the fissioning nucleus for all deformations 
along the fission path. The double-humped barrier (Fig. 1.1') 
which resulted from such calculations has been successful in 
explaining some of the above mentioned characteristics of 
fission. 

2.2 T he Unified Model of Fission 

Before going into the formal discussion of this 
model it is necessary to understand the important feature 
of 'fission channels' , especially for low excitation 
energies of the compound nucleus. In the formal theory of 
nuclear reactions, a channel is defined according to the 
quantum numbers of the asymptotic wave function of the 
system. Thus an entrance (exit) channel is specified by the 
state of the projectiles (products) before (after) any 
interactions occur. Bohr and Wheeler (1939) first intro- 
duced the concept of fission exit channels in terms of saddle- 
point configurations which are energetically available. By 
fission exit channels one means a well-defined type of 
fragmentation of the nucleus into two definite nuclear 
species in given excited states. The number of such 
channels would be extremely large (Wilets, 1964). Therefore 
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the variation in fission width should he practically constant 
from reasonance to resonance as are the radiative capture 
widths. (Michaud on, 1973). However, experiments indicate 
large fluctuations in fission widths. 


The channel theory of A. Bohr (1956), also called the 
Unified model, gives a straight forward explanation for the 
large fluctuations in fission widths. Before discussing 
this theory it would he instructive to discuss the concept of 
'open' fission exit channels as first introduced hy Bohr and. 
Wheeler (1939). The number, H, of such 'open' channel is 
given hy 



where <^)and (j)} are the average fission width and spacing 
of the reasonances, respectively. An effective number 
of exit channels is defined for each spin and parity J 71 as 


^ K eff ^ jit 
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where Ep. is the sum of penetrabilities P. for all the 
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channels i having the same spin and parity J n , and 

-g i 

P ± = [l+exp (-2u - s r-)]' 1 


is the height 


of the fission barrier associated with the 


transition state i. 
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% has been calculated for low energy neutron induced 
235 239 

fission of U and Pu using experimental data on integral 

fission widths. (Blons et al. , 1970, 1971). The values of 

(J^eff) 0 ^^ a ^- ne< ^ ^ or eaG ^- of these resonances is usually less 

than or equal to one. This indicates that the effective number 

of fission exit channels in the case of low energy fission is 

very small. Por example in the case of ^^Pu (l/2 + ) resonances 

induced by s-wave neutrons have spin and parity 0 + and 1 + . The 

value of IT ££ for each of these two families of resonances is 

close to unity (Michaudon, 1973) indicating that in s-wave 
239 

fission of Pu there are two 'open' channels corresponding 
to the two resonances 0 + and 1 + . This is the basic idea behind 
Bohr’s channel thoory. Bohr had, however, proposed his theory, 
outlined below, before these calculations were made. 

Por excitation energies of the compound nucleus not 
too far above the fission threshold, the nucleus, in passing 
over the saddle point, is " cold" since a major part of its 
energy is bound in the potential energy of deformation. The 
quantum states available to the nucleus at the saddle point 
(the fission channels) are then widely separated and represent 
relatively simple types of collective motion of the nucleus. 
These states are expected to resemble the ground state 
excitations of the deformed nucleus, that is, the channels 
at the saddle point form a spectrum similar to the low lying 
states of the stable deformed nucleus in its ground state. 
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Pig. 2.1 shows the spectrum of nuclear levels for 
an even-even nucleus as a function of a symmetric deforma- 
tion parameter (for the sake of simplicity an LDM type of 
barrier has been chosen). It is assumed that the nuclear 
shape remains axially symmetric during the passing of the 
saddle point. The channels can then be characterized by the 
quantum number K , representing the component of the nuclear 
angular momentum J, along the symmetry axis. For even-even 
nuclei the lowest state of the nucleonic structure has K=0, 
corresponding to a paired nucleon configuration. With this 
intrinsic state is associated a rotational band with energies 

E J,K = E K + It j ( j+1 ) 

where J is the total angular momentum (spin) and I is the 
effective moment of inertia. This latter quantity depends 
on deformation. For the ground .states of the very heavy 
elements the value of ft / 21 is about 7 keV and at the saddle 
point the value would be smaller because of larger deformation. 

For nuclei whose shapes possess reflection symmetry 
the spectrum, for K=0, contains only rotational levels with 
even J values % 0, 2, 4 ..... all of which have positive 
parity. The separation of these transition states comprising tb 
the lowest rotational band would/be smaller than observed in 
the heavy nuclei (ground state) spectra - the factor ft / 21 

*|* — J* 

being smaller. The 2 state is about 11 keV and the 4 state 




nucleus 
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about 37 keV above the 0 + threshold (Lynn, 1968). The 
rotational band also contains odd J values having negative 
parity. These states correspond to excitations of the vibra- 
tional states of the deformed nucleus which do not possess 
reflection symmetry [the 'sloshing' mode of vibration, 

Wheeler (1956)]. These states are antisymmetric with respect 
to rotations about 180 deg and hence give rise to negative 
parity states such as 1”, 3 . They are displaced with 
respect to the positive parity states by an amount hw, where 
w is the frequency of 'tunnelling' motion between the mirror 
shapes of the preferred asymmetry. The higher the degree of 
asymmetry the smaller the w. These low lying (l“) excitations 
of the nuclear ground state have been observed to occur in the 
even isotopes of Ra and Th with energies in the region of 200 
to 300 keV (Stephens et al. , 1954). The energies increase as 
one goes to heavier elements. 

At the saddle-point shape one expects even-even 
nuclei to have a lowest state of J =0 and close lying 
collective excitations of 2 + , 4 + ... type as well as states 
like 1 _ , 3”".... at somewhat higher energies. 

Apart from these collective rotational excitations, 
the nucleus possesses states corresponding to the excitation 
of the nucleonic configuration i.e. single particle excita- 
tions. Since these require the breaking of a nucleon pair 
there will be a significant energy gap (^1 MeV) between the 
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lowest configuration (K=0) and the first excited configura- 
tion, After this gap the spacing between intrinsic 
excitations is only of the order of a 100 keV, corresponding 
to the average spacing of individual particle orbits. With 
each such intrinsic excitation, characterized by a definite 
-K, there is associated a rotational band with J=K, K+l, K+2... 
and both parities. 

In odd-A nuclei the lowest K value is given by the 
component of angular momentum of the last odd nucleon in its 
lowest binding state. At saddle point this lowest K value 
differs in general from the nuclear ground state. The 
spacing between the states of the last odd particle is of 
the order of a few hundred keV, and the intrinsic excitations 
in odd-A nuclei involve no energy gap similar to that 
characterizing the paired configuration in even-even nuclei. 
With each .particle configuration is associated a rotational 
band with J = K, K+l.... and both parities. Similarly in 
odd-odd nuclei the intrinsic excitations would be closely 
spaced . 

It must be mentioned that whether the single-humped 
or double-humped barrier is considered for the purpose of the 
Unified model the conclusions are not altered. As shown in 
Pig. 1.1 both saddle points (A and B) have intrinsic channels 
associated with them and hence the channel which would be 
relevant to the fission process would be the one corresponding 



28 


to tlie higher hump. For light actinides the inner harrier 
A is lower than the outer a harrier B whereas the opposite 
situation prevails for heavy actinides. 

2.3 low-energy fission and the characteristics of transition 
st ates 

Slow neutron capture by an even-odd nucleus of spin 
J Q leads to a compound nucleus of even-even type with 
J = J Q + 1/2 and with the same parity (if 1=0) as the 
target nucleus. As explained in the last section the spectrum 
of states at the saddle point (the transition states) will 
contain only one of these spin parity combinations in the 
rotational band associated with the lowest nucleonic confi- 
guration K = 0. Therefore, the fission threshold is expected 
to differ appreciably (~'l MeV) for the two types of compound 
nucleus levels formed. Hence, in slow neutron fission, the 
transition state is essentially well defined and is expected 
to influence some of the characteristics of the process such 
as (i) angular distribution of the fission fragments (ii) 
fission widths of the resonances (iii) mass distribution of 
the fission fragments (iv) kinetic energy of the fragments 
(v) average number ^ of prompt neutrons emitted and (vi) 
emission probability and energy spectrum of light charged 
particles (1CP) emitted in fission. These characteristics 
are discussed in section. 2.5. 
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In the discussion so far only neutrons with 1=0 

(s-wave) neutrons have been considered. However, the 

fission cross-sections for different partial waves differ 

significantly depending upon the energy of the incident 

neutrons. Rae et al (1958) have calculated the neutron 

235 

capture and elastic scattering cross-sections of U 
upto 1.1 MeV on the basis of the statistical theory of 
nuclear reactions. In the course of these calculations, the 
partial wave fission cross-sections, upto 1=3, were computed 
and these are shown in Rig. 2.2. The curves are dashed 
above 500 keV because above that energy they depend on the 
assumptions used concerning the low-lying level structure 
of 255 U„ It is seen that the p-wave (1=1) cross-section 
rises rapidly with energy and becomes equal to the s-wave 
(1=0) cross section for a neutron energy of 100 keV. Rrom 
there on to 500 keV, the p-wave component represents at 
least 50 X of the total fission- cross -section;, the d-wave 
having risen to equal the s-wave component by 500 keV. 

It follows, therefore, that depending upon the 
incident neutron energy different channels would influence the 
fission process to differing extents. Ror example, in s-wave 
fission of 2 ^U(7/2“) the states accessible to the nucleus 
would be 3~" and 4” and when the energy is increased so that 
p-wave interactions become probable , the levels of character 
2 + , 3 + » 4 + and 5 + would also be accessible. The relative 



(barns) 



E n (keV) 


3 ?i<r. 2.2, Calculated partial wave fission 

cross sections for neutron induced 
fission of 235y. [After Cuningh&me 
et al. » 
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predominance of either of these states i.e. odd or even 
parity states, would depend on the relative predominance of 
the partial fission cross-sections in the neutron energy 
region under study. Hence, if the neutron energy spread is 
not too large it should be possible, in principle, to 
selectively excite one or a few of these states and study 
the characteristics of fission resulting thereby. As the 
neutron energy increases the channel density also increases 
(since the level spacing decreases) and thus many channels 
would participate, and the channel effects are expected to 
average out . 

It should be noted that the assumption that the 
fission characteristics are influenced by the properties of 
the transition states depends on the fact that the scissioning 
nucleus retains some information (over and above the rigidly 
conserved quantum numbers) as to which channel it passed 
through. In order to further understand this it is necessary 
to consider the descent of the nucleus from saddle to 
scission. 

2 • 4 Motion from saddle to scission; Adiabatic and Statistical 
Models 

The excitation energy of the compound nucleus may be 


represented as 
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E* = E. + e 
1 

where is the excitation energy of the i channel above the 
ground state and e is the fraction of E stored as the 
potential energy of deformation (Pig. 2.1). As has been 
mentioned, e would constitute a major fraction of E , in 
slow neutron fission, such that the nucleus is 'cold' at the 
saddle point and would be of the order of a few tens of keY. 

Most models of fission can be divided into two groups 
depending on whether or not the motion from saddle to 
scission is assumed to be adiabatic with respect to the 
particle degrees of freedom (Yandenbosch and Huizenga, 1973). 

If the change in deformation during descent is sufficiently 
slow or the coupling to internal degrees of freedom (residual 
interaction) sufficiently weak so that the single particle 
degrees of freedom can easily readjust to each new deformation 
(i.e. no nucleonic excitation) as the distortion proceeds, 
an adiabatic approximation may be valid. In this situation 
the decrease in potential energy from saddle to scission 
remains in the collective degrees of freedom at scission and 
appears primarily as kinetic energy associated with the 
relative motion of the nascent fission fragments. Hence the 
motion will become faster as fission proceeds. If, on the 
other hand, the change in deformation is somewhat faster or the 
coupling somewhat stronger there will be transfer of collective 
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energy into nucleonic excitation in a manner analogous 
to viscous heating. If there is sufficient non-adiabatic 
mixing of the energy among the single particle degrees of 
freedom by the time scission is reached, then a statistical 
model may be a reasonable approximation. The transfer of 
collective energy to the internal degrees of freedom can ■> 
be thought of as arising due a nuclear viscosity. In this 
case less energy will be available for kinetic energy of 
motion of the fragments hence the descent would become 
slower . There is thus a feedback mechanism tending to 
prevent either approximation from becoming increasingly valid 
as the motion from saddle to scission proceeds. 


At the other end of the time scale is the 'sudden' 
(impulse) approximation. Here it is assumed that the descent 
is so rapid that the nucleonic motion cannot readjust at all 
to the change in potential. A condition for this approxima- 
tion to be valid is that the descent time is short compared 

—2 2 

to nucleon times .( ^ 3xl0~* sec.). An estimated descent 
-21 

time of '•'•'3x10 sec has been derived from 1IM calculations 


(Hill, 1958). 

» io). 
n 


(Blyumkina et al. , (1964) had obtained 
Thus the descent is slow enough for the impulse 


approximation to be invalid. However, the actual tearing of 


the neck and the collapse of the nascent fragments may be 
sufficiently rapid for this approximation to be valid 


at scission. 
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It is difficult to make theoretical estimates which 
are sufficiently reliable to decisively indicate whether an 
adiabatic approximation is valid. Wilets (1964) has examined 
this question from the point of a level crossing model 
(Hill and Wheeler, 1953) and concluded that although the 
process is unlikely to be adiabatic, it is also unclear 
whether the non-adiabatic mixing is sufficiently strong to 
assume the validity of the statistical model. Wilets has 
emphasized that the large energy gap associated with pairing 
may cause fission of even-even nuclei at low excitation 
energy to be mostly adiabatic. In reality the descent from 
saddle scission is undoubtedly very complicated, with near 
adiabaticity occuring for some degrees of freedom and 
sufficient non-adiabatic mixing for others. Experiments on 
low energy neutron induced fission might, however, indicate 
the extent to which channel effects show up after scission. 
Such information could help decide which of the descent 
models is most applicable. 

2 . 5 Channel Effects in low-energy neutron induced fission 

2.5.1 A ngular distribut i on of fission fragments 

Perhaps the most striking phenomenon explained by 
the channel theory of Bohr (1956) was the angular distribu- 
tion of fission fragments, using the simple symmetric top 
model and the quantum numbers characterising the fission 
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channels. The angular momentum of the compound nucleus J, 
and its projection on the symmetry axis K are schematically 
shown in Fig. 2.3. The direction of motion of the fission 
fragments is along the symmetry axis. .For relatively small 
neutron energies the angular momentum vector remains 
predominantly perpendicular to the symmetry axis, (as also to 
the incident direction) hence only smaller K values play a 
role and the distribution of fission fragments is peaked in 
the forward direction. This explanation, of course, implies 
that during descent, the value of K does not change. For 
larger incident energies, a larger number of channels may 
enter the picture, nevertheless, a systematic tendency for 
anisotropy in fragment emission results from the fact that 
the values of M (the projection of J on the incident neutron 
direction) and K are not distributed from - J to +J but are 
concentrated towards numerically small values. Such aniso- 
tropies have been experimentally observed (eg. Nesterov et 
al. , 19675 Bekhami et al. , 1968$ Nadkami, 1969) confirming the 
results of the channel theory. Anisotropies in photo-fission 
(Baerg et al. , 1959) and in fast neutron- induced fission 
(Brolley et al., 1954) and other particle induced fission, 
such as (d ,pf ) and (t ,pf) , reactions have been explained 
similarly (Bohr, 1956; Wilets and Chase, 1956; Griffin, i960). 



an^uler momentum coupling scheme for a deforms 
nucleus. The vector' J defines the total anyul 
uioaontun . The quantity M is the- component of 
on the apace-fixed. Z axis (the bean direction) 
The quantity K is the component of J along the 
nuclear symmetry axis* The collective rotatio 
angular momentum R ia perpendicular to the nuc. 
symmetry axis [After Vandenbosch and Huissenga* 
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2-5.2 Fission Widths 

Since slow neutrons absorbed by odd mass nuclei form 

states with either of two spins and a parity well defined by 

the parity of the target, one expects a division of observed 

resonance fission widths into two groups with different 

average values of depending on the spin of the resonance 

state. As an illustration we can examine the well studied 
239 

Pu fission resonances. Resonances induced by s-wave 

P3Q jr -i_ 

neutrons in Pu have J =0 or 1 , The low spin value 
of these resonances leads to relatively large spacings and 

'T r -J- 

_ neutron widths. For the J =0 fission resonances, there 
is at least one fully open 0 + fission exit channel, that of 
the ground state, since it lies about 1 MeV below the 
neutron separation energy, S n . (The K = 0 channel is 
the main fission mode). In contrast, no simple 1 + collective 
excitation seems to exist below S n +1 MeV where two-particle 
states start to be excited. 

The large energy difference between the lowest 0 + 
and 1 + transition states is actually reflected in the fission 
width distribution obtained experimentally (Blons et al. , 
1970). Two families of resonances are clearly visible and 
these two groups of narrow and wide fission resonances have 
been shown to have = 1 + and J 71 = 0 + respectively. This is 
in excellent agreement with the predictions of the channel 
theory. 
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2.5.3 Mass distribution of fission fragments 

A suggestion stemming from the nature of the saddle 
states is that for resonance fission the mass asymmetry of 
the fragments may also depend upon the spin and parity of 
the resonance state. A tentative mechanism for such a 
dependence was suggested by Bohr (1956), Wheeler (1956) 
and discussed by Wilets and Ohase (1956). One assumesthat 
the mass asymmetry may be influenced by the state of 
asymmetric vibration at the barrier through which the 
nucleus passes en route to scission. In particular, fission 
via the first excited . state of the asymmetric vibration 
could show a decreased probability for symmetric mass 
division relative to fission via the ground state of this 
vibrational mode. Hence one would expect that fission 
proceeding via the even parity band e.g. states based on 
0 + , would show a more symmetric mass distribution (higher 
valley in the mass yield curve) than that proceeding via the 
odd parity band e.g. states based on 1”. Odd-spin even- 
parity (1 + , 3 + ....) or even-spin odd-parity (2", 4”....) 
levels, which belong to neither band, would contain wave 
functions of both the symmetric and antisymmetric type (in tho 
vibration parameter) so that their mass distributions might 
be intermediate in character between the other two groups. 

A very large number of radiochemical ond other measurements 
of the peak-to- valley ratio have been carried out for a 
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large number of resonances in low energy (upto 600 eV) 

p rz q 

fission of "^Pu (e.g. Cowan et al. , 1966; Burgus , 1962; 

Melkonian and Mehta, 1965; Ryabov et al. , 1972) and of 

(Regier et al. , I960,* Cowan et al. , 1961; 1963, 1970; 

Simpson et al. , 1971; Ryabov et al. , 1971; Melkonian and 

Mehta, 1965). The results indicate that this ratio, R, 

does change from resonance to' resonance and that they 

appear to fall into two symmetry classes characterized by 

239 

different average values of R. In the case of Pu, the 
average values of R for these two groups were in the ratio 

4. -j- 

3s5 corresponding to the two spin states 0 and 1 . As 

expected from theory, low-R and high-R groups belong to the 

0 + and 1 + spin states, respectively. There was also a 

correlation between the R values and The widths were 

larger, on the average for the low-R group of resonances. 

235 

In the case of U, the distribution in R values was 
narrower but they also appeared to be composed of two groups 
of resonances though the separation was not as clear as for 
259 Pu (Michaudon, 1976). No correlation was observed 
between R and 

In s-wave (1=0) fission the levels accessible in the 
case of 2 '^U(7/2~) are 3” and 4” and for 2 ^Pu(l/2 + ) are 0 + 
and 1 + . In p-wave (1=1) fission the levels are 2 + , 3 + , 4 + 
and 5 + for 2 ^U and 1“ and 2” for 2 "^Pu. Using the 
considerations developed in the channel theory one would 
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expect more symmetric fission (smaller B) for the even 

parity states, in particular the even spin states i.e. 

2 + , 4 + for 235 U and 0 + for 2 ^Pu, than for the odd parity 
states. Such effects have been observed (Cuninghame et al. , 
1961, 1964; Mehta et al., 1967). Cuninghame et al. , (1961, 
1964) measured these ratios for fission by neutrons in the 
energy range between 6 5 keV and 1 MeV with emphasis on the 
region around 125 keV. 

235 

The results on U fission indicated a large 
increase in R below 200 keV incident neutron energy. The 
measured value of B at 125 keV was approximately 45% above 
the thermal neutron value. The increased ratio in the 
energy region of 60 to 125 keV especially at 125 keV is 
quite contrary to what one would expect from channel theory 
considerations. The explanation of these anomalous results 
was not conclusive and some of the suggestions to account for 
thes©.-were ( i) in thermal fission the mass distribution is 
only averaged over a mixture of a few levels of opposite J 
and statistically this type of fission may not be truly 
representative of s-wave fission; as the energy was increased 
the relative number of fissions via even parity levels increa- 
sed hence the' proportion of symmetric fission increased, 
however, the large increase in the 125 keV region remained 
unexplained (ii) the abiabatic model does not hold good 
so that between saddle and scission there is a mixing of 
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adiabatic model states j however, in order to provide lower 
symmetry for p-wave fission, one would have to assume mixing 
between levels of different K. In other words, it would be 
necessary to violate K-conservation which appeared so 
necessary to explain angular distribution data. 

239 

The -Tu results indicate an increase in ratio in 

the region 30 keV to 200 keT above the thermal value, the 

increase being very significant around 30 keV. Beyond this 

region the ratio is more or less the same as the thermal 

value. The region (30 keV to 200 keV) of increase is the 

region where p-wave fission competes with s-wave fission 

and the increasing asymmetry (large R) would be due to the 

relative predominance of negative parity levels populated 

due to the p-wave interaction. These results are quite 

consistent with the predictions of the channel theory unlike 
235 

the " U results. 

2.5.4 Kinetic energy of fiss ion fragments 

The kinetic energy of the fragments arises out of the 
Coulomb repulsion at scission with the possible addition of 
prescission kinetic energy which depends on viscosity 
effects. Measurements of kinetic energy distribution in 
resonance neutron induced fission of (Melkonian and 

Mehta, 1965? Moore and Miller, 1965) and 2 ^%>u (Melkonian 
and Mehta, 1965) showed that it was possible to separate 
the variations in the kinetic energy distribution into two 
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groups corresponding to the two different J values in s-wave 

fission. Similar measurements on the average total fragment 

235 

kinetic energy in U fission in the incident neutron 

energy region upto 2 MeV (Bolshov et al. , 1968$ Blyumkina 

et al. , 1964) have indicated that the average total fragment 

kinetic energy shows definite structure around 200 keV - it 

drops by about 0.6 MeV in this region. This drop in kinetic 

energy has been explained as resulting due to the predominance 

of the p-wave interaction leading to increased excitation of 

the fragments. (The number of prompt neiitrons emitted (V ) 

is a measure of this excitation). The positive parity 

states populated by the p-wave interaction lie about 0.6 MeV 

below the negative parity states populated by the s-wave 

interaction. Measurements on • U (Okolovich and Smirenkin, 

1963$ Boldeman et al. , 1976) which has a ground state parity 

235 

opposite to that of U$ showed an opposite effect in average 
total fragment energy around 200 keV. These results go to 
show that channel effects persist even in gross quantities 
such as average total fragment kinetic energy. 

2 . 5 . 5 Prompt neutron emission 

A major portion of the energy not appearing as 
kinetic energy is dissipated by emission of prompt neutrons. 
Thus the average number of prompt neutrons emitted ( v ) 
would be a direct measure of a fairly large amount of 
energy stored in the nascent fragments at scission. This 
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energy may either he in the form of internal (excitation) 
energy, or be stored in the form of deformation energy and 
converted into excitation energy as the fragments collapse 
to their equilibrium shapes. The linear dependence of HF 
on incident neutron energy was postulated by Fowler [Fowler’s 
hypothesis, Leachman (1956)]. According to this hypothesis v 
increases linearly with neutron energy as a result of 
increasing excitation energy - the average total kinetic 
energy being independent of neutron energy. Andreev (1961) 
had predicted that in the energy region E n <l MeV the fission 
channels would effect the distribution of energy between 
kinetic energy and excitation energy of the fragments leading 
to deviations from Fowler’s hypothesis. This prediction is 
based on the fact that because discrete channels are available, 
the average kinetic energy must increase linearly with 
incident neutron energy and y” remain constant for fission 
through a particular channel. When a new channel opens, the 
kinetic energy decreases in a step and V increases propor- 
tionately. Andreev (1961) had assumed that the energy that 
goes into deformation of the fragments^as the energy is 
increased from one channel to the next (higher barrier) f is 
converted into kinetic energy of the nascent fragments, 
whereas Fowler’s hypothesis requires this deformation energy 
to go into nucleonic excitation. The validity of these two 
lines of thought are, we see once again, connected with the 
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nature of the descent from saddle to scission. 

The variation of V from resonance to resonance has "been 
studied for the fission of 23 ^Pu (Shapiro, 1968 j Weinstein et al., 
1974) and 235 U (Ryabov et al. , 1972; Trochon et al. , 1973). 

These studies showed definite variations in V which could he 
correlated, though not quite consistently, with the different 
values of J n (Michaudon, 1976). However, with the establishment 
of the existence - of the (n, yf ) reaction from results obtained 
on measurements on >' and (the mean fission y-ray energy) 
(Prehaut and Shackleton, 1974)? the variation in 77 could be 
explained without invoking the influence of the transition 
states. The conclusion is that there is a weak spin effect on 
2> but most of the variation in v is due to the effect of the 
(n, yf) reaction and not to the fission channels. 

Measurements on V for higher incident neutron energies 
(upto 1 MeV) in the fission of 233 U (Blyumkina et al. , 1964? 
Meadows and Wahlen, 1967) and 233 U (Blyumkina et al. , 1964; 
Boldeman et al., 1976) have indicated the presence of 
structure around 200-300 keV neutron. energy. The increase 
around 200 keV relative to thermal energies was correlated 
with the decrease in average total kinetic energy of the 
fragments (~ , 0.6 MeV) in this region. The interpretation of 
the energy division was that the energy associated with 
collective motion (eg., rotation) was weakly coupled to the 
nucleon degrees of freedom and hence this energy went mostly 
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into kinetic energy of the fragments. Thus, whenever a new 

channel opens, there would be change in kinetic energy and 

255 

'between' channels v would increase. In the case of U, 

since the even parity states (populated by the p-wave inter- 
action) lie~0.6 MeV lower than the odd parity states, the 
kinetic energy would also register a drop in the region 
where p-wave fission becomes important (-^200-500 keV) and the 
' extra' energy made available would go into fragment excita- 
tion and subsequently result in an increase in v. Beyond 
600 keV channel mixing effects would probably result in the 
averaging out of v and as observed. 

2.5.6 Light charged particle (LCP) emission 

For the scissioning nucleus to be able to emit a 
light charged particle (LCP) it is necessary that a certain 
amount of energy be concentrated on a degree of freedom which 
facilitates LCP emission, for example Halpern (1963) has 
calculated that for alpha particles (LEA) this energy is about 
29 MeV, Since the distribution of energy between the various 
degrees of freedom eg., kinetic energy, fragment excitation 
etc., at the saddle point, is influenced by the properties 
of the saddle states it would be reasonable to expect that 
the emission probability (yield) and average energy of LCP 
may also be affected by these states. 
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Measurements on LRA emission probability (Pjj^) have 
been carried out for fission induced by neutrons in the 
resonance energy region (Deruytter and Ileve de Mevergnies, 

1965; Melkonian and Mehta, 1965; Schr'oder et al., 1965; 

Wagemans and Deruytter, 1972). The results indicate that the 
variation in P^-^ (as indicated by the LRA yield) and average 
energy of IRA is in correspondence with the J value of the 
compound nucleus resonance i.e., fission channel. Measure- 
ments with fast neutrons (Drapchinski et al. , 1964) and with 
fast charged particles (Thomas and Whetstone, 1966 $ Loveland 
et al., 1967) have also been made. At high incident energies 
the probability for second and third chance fission increases, 
hence the observed Pj^ would be an average over several 
nuclear species fissioning at different excitation energies 
leading to uncertainties in Pjj^* Hence, it would be necessary 
to measure at those energies where only first chance 

fission occurs. In the case of and 2 ^Pu this (incident 

neutron) energy range lies below 5 MeV. 

Investigations on P^-^ in neutron induced fission of 
upto 4 MeV, in intervals of 1 MeV (Hadkarni and Kapoor, 
1970) .showed that P^^ was equal to the thermal neutron value 
within the statistical error of 10 Subsequent measure- 

ments with neutrons in the range 0.75 to 1.75 MeV also gave 
the same result (Hadkami et al., 1975) however, the low- 
energy component of the spectra showed variation with neutron 
energy. In the energy range above 500 keV channel mixing 
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effects could result in an average value of close to 

the thermal value. Hence, in order to determine the 
presence of structure in yield or average energy one would 
have to concentrate on the region from thermal to 1 MeV. 

The present measurements on were for neutron 

induced fission of U (7/2“) and 2 ^%>u (l/2 + ) in the 
energy range from thermal to 1 MeV, with a view to determine 
whether or not channel effects show up in the yield and 
average energy of the LBA. 



CHAPTER III 


EXPERIMENT 


3.1 Exp erim en tal s etup 

The present work was done with the 2 MV Van de 
G-raaff accelerator at the Indian Institute of Technology , 
Kanpur. It is a standard H.V.E.C. machine capable of 
providing positive ion beams upto a total current of about 
100 \xA with an energy resolution better than + 5 keV (with- 
out the stabilization system). A schematic diagram of the 
beam transport system is shown in Pig. 3.1. 

The quadrupole lens (Q-jQg) used 1:0 focus the beam, 
was fabricated locally. The pole gap of this doublet is 
large enough to allow a l^-” diameter ■ beam pipe to pass 
through. The stand on which the lens is mounted has 
facilities for easy up-and-down movement of the lens. With 
this lens it was possible to focus the beam to a spot less 
than 5 mm in diameter at the target located about 6 meters 
away from the -lens. 

The neutron counter used to monitor the neutron 
flux was a standard long counter . (Hanson and Mckibben, 1947) 
It consists essentially of a BP^ counter axially embedded 
in a paraffin cylinder. The exact geometry and other 
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details are described in the reference quoted above. This 
type of counter has a flat response to neutrons above 1 Me-7 
The response to thermal neutrons is about 70% of the 
response to 1 MeV neutrons (Allen, I960). The counter was 
placed at a distance of about 3 meters away from the neutron 
target. It served essentially as a monitor of neutron target 
deterioration. 

The neutron target holder assembly, shown in 
Fig. 3.2, consisted of a water cooled collimator and a 
target holder made of copper. With this assembly it was 
possible to use proton currents upto 50 pA on the target 
without significant target deterioration for periods over 
24 hours. Neutrons were produced by the 'Li(p,n) Be and 
T(p,n) He reactions. The lithium targets, of about 20 keV 
thickness to 2 MeV protons, were prepared by vacuum 
evaporation of lithium metal on to copper backings. The 
tritium targets were obtained from the Isotope Division, 

BARC, Bombay. The tritium targets of activities upto 12 Gi 
consisted of titanium deposited on a copper disc in which 
tritium was adsorbed such that the ratio of number of atoms 
of tritium to titanium was lsl. In order to calculate the 
thickness of these targets in keV the following procedure 
was adopted 

i.i ?. . 7 KFtflR 
CENiK.'.i UBRAlff 

Ace. N o. A. 


Consider a 1 Ci target : 
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M _ Aw 

dt - - AN = 

3.7xl0 10 sec" 1 

n = 

3.7x10-^x12.33x3.15x10^ atoms 

= 

1.44xl0 19 atom So 

Hence the 
number of = 

Ti atoms 

1.44xl0 19 


If this titanium is deposited over an area of 1 cm 2 * then the 
thickness of the titanium layer would be 

1.44xl0 19 x48 , / 2 

— —55- — = 1.14 mg/ cm 

6xicr 9 

The protons lose most of their energy in titanium hence 
using the energy loss of protons in titanium in units of 
MeV/mg/cm~ , obtained from tables, the nett loss for any 
target of specified activity may be calculated. For example 
the energy loss of 2 MeY protons in titanium is 94 keV/mg/cm j 
hence in a 4 Ci target the protons would lose about 130 keV 
i.e. the target thickness of a 4 Ci target to 2 MeV protons 
is ~ 130 keV. 

235 

The fissile targets consisted of deposits of U and 
259 Pu on nickel discs. . The 2 ^U target of 93'. 9 °/ B enrichment , 

obtained from the Isotope Division, BAEC, had a thickness 

2 2 239 

of - 5 mg/ cm^ and an active area of ~4 cm . The y Pu 

targets of 96.97% enrichment , obtained f„rom the Oak Eidge 

2 

national Lab., USA, had a thickness of'- 880 pg/cm with an 

2 

active area of ~ 3 cm . 
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and 410 amplifiers and fed to timing single channel analysers 
(Ortec 420) to get the timing signals. The discriminator on 
the fission side was set to cut off natural alphas ("6 MeV) and 
on the 1CP side the cut off was ~6 MeV in the runs and 

<~'8 MeV in the ^^Pu runs i.e. the cut off was just above the 
highest alpha energy of the calibration source. The timing 
signals were fed to a time-to-puls e-height converter (Ortec 
457, 2t = 0.5 usee.) in the runs and to a slow coinci- 

dence unit (Ortec 409, 2 t = 1 psec.) in the Fu. runs. The 
1CP detector pulses after suitable delay (Ortec 427) were 
gated with these outputs and recorded on a Packard 400 
channel analyser. The fission and alpha count rates were 
simultaneously monitored. 

5 . 2 Procedure 

Neutrons of five different energies in the range 100 to 
1100 keV were produced. The li(p,n) reaction was used to 
produce neutrons upto 500 keV (Gibbons and Newson, I960) and 
the l(p,n) reaction for neutrons between 500 and 1100 keV 
(Coppola and Knitter, 1967). The neutron energy spread was 
calculated using energy loss and angular spread considerations 
(the fissile target intercepted the neutron ’beam’ in a 
50° cone). Thermal neutrons were produced by interposing, 
a 5 cm thick paraffin block between the neutron target and 


the fission chamber. 
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The average proton current used was 30uA. The 
23 ^ 3 

fission rate in the U runs was~2.5xlCr per sec and in 
239 2 

the Pu runs it was ~ 2. 6x10 per sec. The neutron flux 
at the fissile target estimated from these rates was — 10^n/ 
cm / sec. The average 1CP count rate was 1-2 per min. and each 
spectrum was collected for about 20 hours. 

After each coincidence run an ungated (singles) 
spectrum of the 1CP was recorded for a fixed time. This was 
used in the correction of -the coincidence spectra for chance 
coincidences. Each *fast' run was preceeded and followed by 
a thermal run. These thermal runs served as a monitor for 
possible drifts and also provided a normalization for the 
’fast’ runs. The stability of the system was monitored 
frequently using an on-line precision pulser. The bias on the 
detector changes with increase in leakage current because of 
the increase in voltage drop across the series resistance in 
the preamplifier. When the bias drops so low that the deple- 
tion depth is not sufficient to stop the high energy ICPs, 
the spectrum gets distorted. Hence the leakage current in 
the detector was continuously monitored and the bias was 
adjusted when this current increased(due to neutron damage) so 
that the high energy tail of the LCP spectrum was not lost. 
Detector deterioration due to radiation damage was quite severe 
and during the course of the experiment several detectors were 
used. Several coincidence runs (with similar calibrations) 
were added to give the final coincidence spectra. 
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5.3 Data An aly sis 

The LCD spectra were first corrected for chance coinci- 
dences using the singles LCD spectra. The chance coincidence 
rate is given by 


w ? *1 N 2 -1 

N ch ~ 2t T-j_ T 2 sec 

N f -1 

= 2t -f-f sec 

Q 

where 1ST is the number of events in the singles run. 
c 

T is the time period of the coincidence run during 

which the fission rate was monitored, 
s 

T is the time period of the singles run. 


The number of chance coincidences in time T c is therefore 




ch 


= 2t 


a f 


T 


s 


The number of chance coincidences per fission event is 


IT 


H ch - 2 * 


SL 

W 


hence the number of chance coincidences during the coincidence 
run is 


"oh = ^ 


N s u2 

a f 

T S 


or 


n 

§ N ch f i = 


2r m 


,8 
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N 8 . 
i=l a 
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Therefore, if F s . is the number of counts in a parti- 

CC y 1 

cular channel i then ST^ ^ gives the chance coincidence 
correction for that channel. In this manner a channel by 
channel correction was made. The maximum correction (^5%) was 
in the low energy region. 

These spectra were then corrected for energy loss in 
the source, gas and aluminium foil. This was done by 
estimating the average thickness of source, gas and aluminium 
foil encountered by the particles by comparing the observed 
LOP peak energy in the thermal runs with the known value of 
LRA. energy in thermal fission of 233 U (Fadkarni and Kapoor, 
1970) and in 239 Pu (Krogulski et al. , 1969). The shift in 
peak energy represents the average energy loss (-&B) of the 
LEA, with initial energy (E ) corresponding to the known peak 
value, in the source, gas and foil. The average thickness of 
argon corresponding to this energy loss was obtained from 
tables (Forthcliffe and Schilling, 1970) „ Using this value a 
plot of initial energy (E ) versus residual energy (E Q -/iE) 
was made using energy loss tables. It was thus possible to 
determine the energy loss for any initial LRA. energy. The 
energy loss values obtained this way were also cross-checked 
by calculating the energy loss in the source, gas and foil 
step by step and the results indicated that the method 
adopted above was very satisfactory. 
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The experimentally measured energies, which are the 
residual energies, were thus converted into initial energi eS ' 
The final spectra were obtained by adding the counts in 
channels falling within equal initial energy intervals of 
0. 5 MeV each. 



CHAPTER IY 


RESULTS AND DISCUSSION 


4.1 fission 

The LCP energy spectra in thermal and fast neutron runs 
are shown in Pig. 4.1 (Krishnarajulu et al. , 1977). The least 
squares Gaussian fits to the spectra are shown as continuous 
curves. The lower limit of the Gaussian fitting function 
was fixed at 12 MeV since the yield in the region beyond 
12 MeV is predominantly due to alpha particles (see Table 1.1). 
The spectra show more low-energy particles than that given by 
the Gaussian. 

Prom these fits, the most probable value (E a ) and 

standard deviation (a^ ) were obtained as a function of 

a 

neutron energy and these values, tabulated in Table 4.1, are 
shown in Pig. 4.2. It can be seen that, there is some struc- 
ture (•'•'*0.6 MeV) in the neutron energy range of 200 to 700 keV. 
There is no indication of a systematic trend with neutron 

energy. The width, Og , seems to be insensitive to variation 

a 

in incident neutron energy. 

The results on 1CP yield above 6.5 MeV (Table 4.1) 
are .shown in Pig. 4.3(a). There is a marked increase of 
about 20/^ in the range 200 to 500 keV as compared to the 
thermal neutron yield. At higher neutron energies the 
yield tends towards the thermal value. This yield has 



T able 4.1 

2‘55 

Average IRA energies, widths and Normalized Yields ( U fission) 
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Fig, 4,5* "Variation of (a) total LCP yiela above ,5 H 
g (b) alpha yield above 12 MeV and (c) LCP yiej 

between 6.5’ and 12 MeV » with neutron energy. 

^235y fi33i 0n ) ■, 
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contributions from both LRA and other light particles, the 
majority of which would be tritons (Table 1.1). Assuming 
that above 12 MeV the contribution to LCP yield is only due 
to IRA, an attempt has been made to estimate the yields of 
LRA and 'tritons' separately. 

The observed LRA yield for 12 MeV as a function 

of neutron energy is shown in Pig. 4.3 (b). There is an 
increase of about 20 % in the energy region 200 to 500 keV 
similar to the LCP yield results. At higher neutron energies 
the yield is the same as that for thermal energies. Beyond 
1 MeV this yield has been measured (Fadkarni and Kapoor, 1970;' 
and found to be equal to the thermal value within the limits 
of error. 

The yield in the energy region 6. 5 k. PpQp<l 12 MeV has 
contributions from both LRA and 'tritons'. To estimate the 
yield of 'tritons' alone, the low-energy tails of the 
Gaussian fits to the spectra have been subtracted from the 
observed yields to give the 'triton' yield. The extracted : 

' triton' yield for different neutron energies is shown in Pig 
4.3(c). For neutron energies around 200 keV this yield is 
equal to the thermal value but it shows a significant increas 
(2 - 2.5 times) over the thermal value, around 500 keV. ; 

Beyond 500 keV the yield decreases but is always higher than ! 


the thermal value 
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25Q 

4 * ^ __ J Pu fission 

The LCP energy spectra in thermal and fast fission are 

shown in Pig. 4.4. The continuous curves represent the 

Gaussian fits. Since the out-off in this case was larger 
235 

than in the U runs (see section 3.2) the lower limit of 
the spectra are shifted towards the high energy side. The 
fits were obtained in "lie same manner as before, setting the 
lower limit of the fitting function at 12 MeV. 

The most probable energy (P a ) and the standard 

deviation ( 0 - ) as a function of neutron energy (Table 4.2) 
a 236 

are shown in Pig. 4.5. As in U fission, the average 
energy and width do not vary significantly with neutron 
energy and the values are very close to the thermal value. 
There is, however, an indication of a drop in (^0.7 MeV) 
in the region 250 to 300 keV. 

The results on LCP and alpha yield (T^bJe 4.2) shown 
in Pig. 4.6(a), indicate that the variation of yield in this 
case is not similar to or as pronounced as in the case of 
fission. The alpha yield increases slightly (~7%) 
around 250 keV, 600 keV and 1 MeV. One point to note is 
that the alpha yield shows a. drop fi4%) in the range 250 to , 
300 keV. In the other regions of . neutron energy the yield 
is equal to the thermal yield. . 

Since the low-energy cut-off was large (~8 MeV), , 

the low-energy yield data could not be obtained and a part ; 
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of the ’triton’ spectrum was lost. Hence the results on 
'triton' yield have a large statistical error as indicated. 

[Pig. 4.6 (c)]. The results indicate that, within the limits 
of error, the yield upto 500 keV is almost equal to the ther- 
mal value. However, one observes that there is an indication 
that the average yield in the region 200 to 550 keV is larger 

than the thermal value. Beyond 500 keV the yield drops 

255 

below the thermal value unlike the U results where the 
'triton' yield was above the thermal value in the entire 
region beyond 200 keV. This decrease in 'triton' yield is 
not apparent in the total LCP yield because the 'triton' 
emission probability is an order of magnitude lower than that 
of LRA. 

4 • 3 Discussion 

4.5.1 IRA, yield 

25 5 

The increase in LRA yield for C fission in the 

neutron energy region 200 to 500 keV cannot be attributed to 

any excitation energy dependence of LRA yield since the 
« 

change in compound nucleus excitation energy is very small 
and the observed yield does not vary monotonically with 
neutron energy, moreover, beyond 1 MeV previous measurements 
(Hadkarni and Kapoor, 1970) have shown that the yield levels off 
to the thermal value. Thus the increase in LRA yield around 
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200 keV may “be associated with the increase in relative number 
of fissions proceeding via the even parity states populated by 
the p-wave interaction at these energies- Since the positive 
parity states lie lower (^/0.6 MeV) than the negative parity 
states, the extra energy available at the saddle point in 
p-wave fission can go into the degrees .of freedom such as 
deformation, kinetic energy of relative motion of the fragments, 
fragment excitation or into 1CP emission. In binary fission 
Blyumkina et al. (1964) and others have observed a drop in 
average kinetic energy of the fragments and an increase in the 
number of prompt neutrons at these neutron energies as compared 
to thermal energies. The higher LEA yield observed in this 
energy region indicates that the extra energy may also be 
available for IRA. emission. At higher energies higher 1-waves 
contribute to the fission cross-section and this may result in 
the averaging out of channel effects resulting in an IRA yield 
equal to the thermal value. 

Rluss et al. (1972) had measured the triton and LRA 
yields in neutron induced fission of in the range 

thermal to 700 keV. They used both radiochemical and 
particle identifier methods. Their results indicate a trend of 
decreasing LRA yield with increasing neutron energy, however, 
the yield in the entire region was close to the thermal value 
of 2x10”^ per fission. No- structure was observed in this 
region. They had remarked that the reproducibility of 
their results was poorer than the counting uncertainty (~5%) • 
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They had also claimed that their data on indicated a 

u a 

broadening of the spectra at 380 keV, however, an examination 
of their spectra does not indicate this. 

It is not possible to explain the structure in the 

23Q 

Pu' LRA yield in as clear terms especially the small 
increases in the 200 keV, 800 keV and 1 MeV regions. However, 
the drop in yield between 250 and 300 keV may be explained in 
terms of the domination of the p-wave interaction in this 
region. Fission induced by p-wave neutrons would take place 
via negative parity states eg. 1“ state, which lies higher 
(Lynn, 1968) than the 0 + ground state. This means that the 
energy which was available for LCP emission is 'reduced' 
and this is reflected in the drop in yield as also in the 
drop in average alpha energy in this energy region. In the 
region beyond 350 keV various partial waves would contribute 
to the fission cross section and the variation in yield 
possibly reflects the result of the admixture of states 
through which fission takes place. 

It Is therefore clear that since there is no 
monotonic increase in yield with neutron energy, the 
structure in yield cannot be due to a dependence on compound 
nucleus excitation energy. The observed structure in v and 
Ej, in this region for U (eg. Blyumkina et al., 1964 ) has 
been explained in terms of channel effects and there is 
indication that the IRA. yield also reflects these channel 
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effects. More accurate measurements with, better neutron 
energy resolution would have to be carried out in order to 
establish this result more firmly. However, decreasing the 
thickness of the neutron target or reducing the solid anglej 
in order to improve neutron energy resolution, would reduce 
the count rates considerably and each spectrum which was 
collected over a few days would perhaps take weeks. 

4.5.2 1 Triton ' yield 

The method used to extract 'triton' yields is far from 

satisfactory and since the low-ener.gy tails of the spectra 

239 

could not be obtained, especially in the case of Pu, there 
are large errors in the estimates of this yield. However, 
some points are worth noting. The first is that the varia- 
tion in 'triton' yield is quite different from that of LEA 
yield indicating that the emission mechanism may be different 
in the two cases. The second point is that this yield shows 

a distinct increase, outside the limits of the estimated 

25 5 

errors, around 500 keV in the case of U, Whether an 

239 

increase exists for Pu around 300 keV or not would have to 
be established by further measurements (in which the out-off 
is lower). In the 500 keV region both p and d-waves 
contribute to significant extents hence no clear picture 
based on channel effects is possible. 
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The results of Pluss et al. (1972) on triton yield 
indicate that the yield for the entire region from 170 to 
700 keV is about 2 to 2.5 times the thermal value. However, 
their data indicate some structure at 350 keV (particle 
identifier method shows a decrease of ■~20%) and at 700 keV 
(radiochemical method shows an increase of ~50%,) about the 
straight line which can be drawn through their data points. 

The two methods they used, radiochemical and particle identifier, 
do not complement each other in the entire neutron energy 
range and where they do overlap the results (specifically the 
decrease at 350 keV ) are not consistent. They had remarked 
that their data were not sufficiently precise to establish 
the energy dependence of the yields they observed. 

In order to establish the structure in 'triton' yield 
more accurate measurements using AE-E systems for particle 
identification are being carried out. 

235 239 

4.3.5 Comparison of U and Pu results 

The results on the variation of average energies and 
widths indicate that there is no significant difference 
between the two cases. The average energies are equal to the 
thermal values within error except in the region around 
200 .to 300 keV where there is indication of some structure. 

The total 1CP yields indicate some difference 
especially in the region 200 to 500 keV but since the 'triton' 
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239 

yields in Pu have large errors, any comparison of the 
total 1CP yield or 'triton' yield may not he justified. 

However, the LRA yields in the two cases show clear differe- 
nces. The LRA yield beyond 500 keV in levels off to the- 

thermal value , whereas in Pu it is on the average more than 

the thermal value. The yield indicates a fairly 

239 

smooth variation from 200 to 600 keV whereas in Pu the 
yield shows structure ir-1%) in this region. In the region 
below 200 keV the yield variations are identical. Cuninghame 
et al. (1961, 1964 ) had expected the mass asymmetry in binary 
fission to show opposite behaviour in these two cases (below 
200 keV ) but their experiments indicated an essential 
similarity in the two cases. 

The difference in structure in the entire region 

beyond 200 keV cannot be correlated with the parities of these 

nuclei, that is, since their parities are different one would 

expect opposite effects in LRA. yield. It is observed, 

239 

however, that in the region 200 to 350 IroV the Pu LRA 

yield indicates a drop to below thermal values which is 

23 S- 

opposite to that observed in U. In this region therefore, 
the opposite parities of these two nuclei seem to manifest 
themselves. In the region beyond 350 keV the specific level 
schemes of these two nuclei and the effect of higher 1-waves 
contributing to the fission cross-section could possibly 
account for the difference in variation of LRA yield. 



CHAPTER Y 


TRAJECTORY CALCULATIONS 


5.1 Introduction 

The method of trajectory calculations has been 
extensively used to study of the scission configuration. 

This method consists of following the motion of the frag- 
ments (the trajectories) after scission. During this period 
the fragments accelerate under the influence of their 
mutual Coulomb fields and the initial conditions of their 
motion would be given by the state of the nucleus at 
scission eg. the extent of stretching (hence the initial 
separation), initial kinetic energy etc. The final condi- 
tions of this motion would correspond to experimentally 
measurable quantities such as kinetic energies, angular 
distributions etc. Therefore, if the connection between 
the initial conditions and the final quantities was 
established, it would be possible to obtain quantitative 
information on the state of the fissioning system at scission. 

The calculation of these trajectories has been 
carried out classically using a limited but representative 
choice of initial conditions , and the main features of the 
final conditions of motion, such as average energies and 
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angles have been reproduced. The method used in such 
calculations is as follows 2 a set of conditions are chosen 
to describe the scission configuration and using these as 
the initial conditions of motion the trajectories are 
calculated. The 'final* configuration is the one obtained 
after the lapse of a sufficient period of time at the end 
of which the fragments are well separated and the Coulomb 
interaction is consequently very small (^l%of the 
original value). These ’final' conditions are compared with 
experimental results and the intial conditions are modified 
until a set of initial conditions are obtained which produce 
'final' conditions which are in good agreement with 
experimental results. This procedure has a limitation - the 
derivation of such a set of initial conditions only indicates 
that at a given moment the kinematical variables of the 
system coincide with these initial conditions; however, they 
may not coincide with the values of the kinematical variables 
at th e moment of scission . In order to obtain information on 
the moment of scission, the 'earliest' set of initial condi- 
tions which produce satisfactory final conditions should be 
obtained. It is apriori not obvious that an 'earliest* set 
of initial conditions can be found (except for the trivial 
one corresponding to zero kinetic energy) from an examination 
of the experimental distributions. This ambiguity in 
localizing the 'moment of scission' in time is an inherent 
feature in such calculations. It is therefore not 
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surprising that such calculations when carried out for 
binary fission produced a large number of sets of initial 
conditions, all of which could produce satisfactory final 
conditions. However, in the case of ternary fission (eg. 
IRA fission) this ambiguity can be reduced because of the 
nature of ternary alpha emission. The angular distribution 
of the alpha particles is very sharply peaked perpendicular 
to the fission fragment direction (Titterton, 1951) 
indicating that the alpha particle is emitted from the 
region between the fragments when they are still close to- 
gether, at the moment of scission (Tsien, 1948). [A very 
small fraction of alpha particles are emitted from the 
fragments diametrically opposite to those emitted from the 
neck (Piasecki et al. , 1970). It has however, not been 
established that these 'polar' alphas are emitted at the 
instant of scission. ] . This rather accurate localization 
of the moment and point of emission of the alpha particle 
would help fix the moment of scission with less ambiguity. 
It is for this reason that the alpha particle represents a 
powerful tool for the study of the scission configuration. 

A large number of experiments have been carried out 
in which certain features observed in noimal binary fission 
have been compared with the same features in IRA fission. 
The main conclusion reached from such comparisons (eg. 
Fraenkel, 1967) is that, IRA fission has all the essential 
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characteristics of binary fission except for alpha particle 
emission. It is therefore believed that the scission confi- 
guration as obtained from a study of LEA fission would 
represent the general characteristics of the scission 
configuration for normal binary fission too. 

The assumptions made in such calculations are that 

(i) the scission configuration defines the initial conditions 

of the three-body motion, that is, a fissioning nucleus 
is practically simultaneously separated into three 
bodies and the state of motion of the whole system, 
induced by i:he preceeding nuclear interactions, develops 
under the influence of the Coulomb force only. 

(ii) non-relativistic classical mechanics and electrostatics 
provide the necessary theoretical basis for the 
calculations . 

The motion of the three-particles cannot be calculated in 
closed form hence the trajectories are computed numerically. 

5.2 Approximations 

For such a three-body calculation there are 18 
variables, for each of the three particles there are three 
spatial variables describing the position of each particle 
and three conjugate momenta. The condition of conservation 
of linear and angular momenta reduce the number of Indepen- 
dent variables to 12. Since the fragments are confined 
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to the fission axis the number of independent variables 
reduces to 10. 

In order to simplify the calculation, the number of 
dynamical variables used to describe the system must be 
small. The minimum number of variables (as described above) 
is obtained when no internal degrees of freedom are allowed. 
The three particles may be assumed to behave like point 
charges - the point-charge approximation. Strictly speaking, 
the particles should not be taken as points because they 
have nonspherical spatial extension over which their charge 
is distributed. The effect of the charge distribution ought 
not to be neglected especially at the beginning of motion 
when the fragments are highly deformed. These charge 
distributions change with time as the deformation changes. 

An attempt has been made to study the effects of deformation, 
this will be discussed later. 

The point-charge approximation effectively anulls 
the effect of higher multipoles on the fragment fields; 
hence it would be necessary to estimate the error involved 
in using this approximation. For this purpose the fragments 
are taken to be prolate ellipsoids and a calculation of the 
effect of the dipole and quadrupole interactions in comparison 
with the monopole interaction is used to estimate the error 
involved. The potential due to a prolate ellipsoid of charge 
Ze at (r,0) is given by 
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V ( r , 9 ) 


V ^ a2 “ b2 ) n P 2n (Cos 9) 
7g(2n+l) (2n+3 )r 2n+ ^ 


where only symmetric deformations are considered and ’a' and 
'To* denote the major and minor axes of the ellipsoid. If 
is the radius of an undeformed sphere of the same volume, 
then 


2 3 

mat) = xlt 

o 

and the eccentricity can be defined as 


2 

£ 


a 2 -b 


R 


2 


thus 

v(r e) = V e 2n P 2a (Cos Q) 

E o n 4^(2n+l)(2n+3)(r/E 0 ) 2n+1 


considering terms till n=l i.e. only till the quadrupole term 


V(r,0) = 


3Ze 

R 


£ 2 P 2 (Cos Q) 


| ^ L * fc || , ^ | JN ,, ^ 7 T J - . v -^ 

3(r/R o ) 15 (r/R o ) 3 


The first term is the monopole term 


V 


mono 


= ~ and E(r) 


Ze 

7 


This field induces a dipole moment in the other fragment, 
that is, it polarizes the other fragment. The force on 
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each charge is therefore 




Since each charge is hound under the action of this restoring 
force 


F 



mw 


2 


x 


o 


where m is the mass of the charge and w q is the frequency of 
oscillation about equilibrium. Hence the induced dipole 
moment is 


p = ex 


Ze^ 

2 "2 

mw r 
o 



For a fragment with Z charges 


w 



r 2 xM 


where M is the reduced mass of the oscillating system and is 
equal to one fourth the fragment mass. From the location of 
giant dipole resonances in those nuclei which have such 
masses it is known that (Halpern, 1963) 


-hw «s 


2 2 

223 _ Z^e 

h 2 /M r 2 x 


15 MeV 
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with Tt- — 4/3 MeVf 2 j « 160 MeV, r «20 f 


x = 1/20 f 


2 

The energy stored in this polarization is 1/ 2kx which is 
about 0.2 MeV. Since the initial potential energies for 
r = 20 f are of the order of 100 MeV the contribution due to 
the induced dipole moment is less than 1^ of the monopole 
potential and can hence be neglected. 


The second term corresponds to the quadrupole 
moment of the fragment. Clearly, for undeformed spheres e=0 
and this term does not contribute. However for e^0 we have 


V 


quad 


Ze e‘ 

5 R 


P 2 (Cos ©) 
<r/R 0 ) 5 


Ze 

5 


(a‘ 


-b 2 ) 


P 2 (Cos Q) 


The quadrupole moment is given by 

Q = Ze/5 (a 2 -b 2 ) 


At a distance 'r' from the centre of the ellipsoid along 
the axis, the ratio of the quadrupole to monopole potential 
is 


V 


V 


quad 


= 1/5 


/ 2 -.2 \ 
(a -b ) 


mono 


thus at the centre of the adjacent fragment, r=2a, and 
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Y 


Y 


quad 


1 


mono 


4a' 


for e*»0.7, (b/a)^2s; 0.6 for a 


= ^ [l-(b 2 /a 2 )J 
lOf and R»8.5 f 


/. 0.02 

mono 

Thus the interaction energy between the fragments ms 4 % 
greater than it would be if the quadrupole moments were 
ignored. The alpha particle is released midway between the 
fragments i.e. r»a, hence - j - s ^ our "k^ 11183 ^■® ae 

value at the centre of the fragments. Hence the potential 
the alpha particle sees is about 1.1 times that seen without 
quadrupole effects. Since the alpha particle moves away 
quickly it would see an effective time average of the 
quadrupole moment close to the original value. In view of 
the smallness of the induced electric dipole and quadrupole 
moments the effect of higher multipolec can. be neglected. 

These results are similar to those obtained by Halpern (1963). 
To keep the calculation of trajectories simple it is there- 
fore justifiable to consider the fragments and alpha particle 
as point sources of their Coulomb fields. It must be pointed 
out that even in the point-chargp approximation the fragments 
are considered to have a size, in the sense that if the alpha 
particle, in the course of its motion, happened to come 
within a distance of -^8f from the fragment centre it was 


considered to be absorbed. 
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The second major approximation is that of restricting 
the motion of the three particles to a plane, that is, the 
calculation is two-dimensional. The justification for this 
assumption lies in the fact that most experimental measure- 
ments on alpha emission are made in a plane, that is, the 
fragment and alpha particle detectors are normally placed 
in a plane. Hence a comparison of experimental results with 
the results of trajectory calculations using the two-dimensional 
approximation would be meaningful. Halpern (1963) had carried 
out calculations of alpha particle trajectories in three 
dimensions and the results, for. example, on energy and 
angular distribution of alpha particles, compare very well 
with a large number of subsequent calculations (eg. Boneh 
et al. , 1967) using the two-dimensional approximation. Thus 
the restriction to two-dimensions would be a justifiable 
approximation. Using this approximation the number of 
independent variables necessary to describe the three-body 
motion reduces to seven. 

5.3 The motion of the fragments and alpha particle 

The motion of the three particles under the 
influence of their mutual Coulomb fields are calculated 
numerically. The equations of motion are (Boneh et al. * 1967) 

dX . . dU. . 

-a¥ = and m i vf = p ij 
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where X. . is the 3 th coordinate if the i th particle,. U. - 
is the 3 component of velocity and fA is the 3 component 
of the force acting on the i th particle of mass m^. 

These equations are replaced by difference equations ; 


v n+l 

= x n . 

+ 


13 


u n+1 

= u n . 


3-3 

13 


u 11 . 

= u n . 

+ 

13 

13 



U n . .At 
13 


1 P n+1 At 

m i ij 2 

1 -n,n At 

m ± 13 2 


is the average or intern ediate velocity and is used to 

improve the accuracy of the calculation. F^ is the 3 

compoment of the force acting on the i th particle at position 
ttn 

hi - 



The subscript k refers to the other two particles, and the 
superscript n refers to the value of the parameter after the 
n th time interval. 

The time interval is chosen such that when the particles 
are close together, that is, just after scission and the forces 
are rapidly changing; it is small and when the particles are 
far apart their directions etc. , are slowly changing; the 
intervals are relatively larger. Such a procedure ensures 
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that uniform accuracy is maintained at all stages of the 
dynamical process succeeding scission. The total time t 
after n intervals is given by 


hence the size of the n th interval is 


t = t -t , 
n n n-1 


t ( 8 m -e (n - 1)a ) 
o v 7 


= t . ( e a -l) 
n-1 v 


Thus f t ' determines the time scale and 'a 1 determines the 
o 

rate of exponential growth. Roughly ’ t ' determines the 
accuracy at the beginning of the trajectory (t=0) and 'a' 
determines the accuracy towards the end (t=eo). The para- 
meters 1 t Q ' and ’a' are so chosen that the accuracy remains 
more or less constant during all intervals. Tor the present 
calculations the choice was 


t 

o 



sec 


and a 


0.1 


The calculation was terminated after 81 time intervals, that 

is, after 3.3x10”'*'^ sec by which time the distance between 

2 

the particles would be ~10 times their initial distance and 

-2 

hence the potential —10 times the initial value. The final 
energies are given by 




E i (t n } .= 4-1/2 m i [v ij (t n )]' 

J 
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.and the angle the alpha particle makes with the fission axis 

“i rn 

, v a (t ) 
o (t ) = tan 1 “ 

“ " < <*n> 


The calculation rigorously conserves linear and angular 
momentum of the system. The maximum percentage error in 
value of total energy after 81 time intervals was estimated 
to be'"' .15 % . 


5*4 Choice of variables 

The seven initial dynamical variables used to describe 
the scission configuration are shown schematically in Pig. 5.1. 
They are ; 

(i) The mass ratio R = m /m- r which is approximately equal 

XI -iJ 

to the charge ratio. However, the charges were 
determined using the relations (Mukherji, 1969) 


Z H = m h/ 2 *587 and z j J = z tot~ Z H 


Z^ 0 j. is the charge of the nucleus minus the charge of the 

alpha particle. The present calculations are for the nucleus 

252 , 
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Of 


and thus .j. = 96 


(ii) The line joining the two fragments is defined as the 

fission axis (x-axis) and the origin is the zero field 
point. The x coordinates of the fragments are inversely 
proportional to their charges except for the correction 
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introduced due to the alpha particle. The coordinates 
are denoted by D, the separation between the charge 
centres of the fragments. 


(iii) The x coordinate of the alpha particle X q 

(iv) The y coordinate of the alpha particle Y q , which 

indicates that the alvha. particle is created off the 

fission axis at a distance Y . 

o 

(v) The initial velocities of the fragments in they 
direction are zero since their motion is confined to 
the fission axis and the velocities along the x axis , 
v£ and v° are combined to give the intial total frag- 
ment kinetic energy, E^. 

2 2 

Ep = 1/ 2 + ni Xi V Xj ) 


(vi) The initial energy of the alpha particle, E°. 


E° = 1/2 m v° 

a ' ■ a a 


(vii) The initial angle between the direction of motion of 
the alpha particle and the direction of motion of the 
light fragment, 0?. 

In most calculations the dhoice of dynamical variables is the 
same as the one used here. However, some calculations have 
used t, the time of emission of the alpha particle, as a 
variable. This is equivalent to a different choice of initial 
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velocities for the fragments or alternatively a different 
separation between them. 

3 • 5 Review of prev i ous calculations 

The calculations that have been made so far can be 
divided into two categories; the ones whose initial conditions 
lend support to the predictions of the statistical theory of 
fission (Pong, 1969) and those whose intial conditions 
support the predictions of dynamical theories such as the 
adiabatic model or sudden snap model (Halpern, 19 63 , 1965, 

1971 ). Broadly speaking the statistical model demands that 
the initial energies be low (~1 MeV) for the fragments and 
~0.5 MeV for the alpha particle) whereas the adiabatic model 
predicts large values for the initial energies (~30 to 40 
MeV for the fragments and ~3 MeV for the alpha particle). 

The former condition implies that at scission the particles 
are almost at rest whereas the latter implies that at 
scission the particles are already moving with a considerable 
fraction (~20%) of the final velocities. One general 
similarity in a large majority of these calculations is 
that most of the initial variables had fixed starting values and 
the final values (distributions) were generated by varying 
one or two of these variables. These calculations are 
briefly reviewed below (see also Peather, 1969; Asghar et al. , 
1971; Vandenbosch and Huizenga, 1973)- 
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5.5.1 Calcul ations supporting the Statistical model 

According to the statistical theory, the fission 
fragments have about 0.5 MeV initial kinetic energy at the 
moment ox scission associated with their three degrees of 
freedom. Hence for the alpha particle a value of 0.5 MeV 
is chosen as its initial kinetic energy (Ertel, 1968). An 
extreme case of this model would be one where the initial 
kinetic energies are zero. 

Tsien (1948) assumed the initial kinetic energies 
of the fragments and alpha particle to be zero and 
discussed their initial coordinates in terms of two para- 
meters; the distance of the alpha particle from the fission 
axis and parameters defining the change of radii of the 
fragments as functions of excitation energy. The velocity 

distribution calculated for the alpha particles emitted 

?56 

in ternary fission of U was in agreement with experiment 
in the high energy region. 

G-eilikman and Khlebnikov (1965) had performed these 
calculations using the droplet model. Nuclear deformations 
upto the term were -considered. The initial alpha 
energy lay .in the range 0 to 2 MeV. The best agreement 
between the calculated and experimental results, on energy 
and angular distribution, was observed at an energy E°=l,o 
MeV. The significant point about the method used was the 
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absence of the point-charge approximation. The deformation, 
however, was assumed to be independent of time, that is, 
the change in ag and during motion was neglected. The 
mean values of the final energy and angle agreed very well 
with experimental results, however, the calculated widths 
were smaller due to the rapid fall in the distribution at 
about 18 MeV. 

Ertel (1968) carried out trajectory calculations 
using initial energy values of 1 MeV and 0.5 HeV for the 
fragments and alpha particle, respectively. He was able to 
successfully reproduce the experimental results (Eraenkel, 

1967) on the angular correlation of the alpha particle, 
that is, the variation of 9 -q with R. 

Vitta (1970, 1971) used initial energies as predicted by 
by the statistical theory and was able to satisfactorily 
reproduce the final kinetic energy distribution. The final 
angular distribution could be reproduced if delayed emission 
was assumed. He also obtained the variation of the final 
alpha kinetic energy and angle with T as well as its (E ) 
variation with initial kinetic energy. 

Eong (1970) extended Ertel' s calculations and 
showed that the angular correlation of the alpha particle 
obtained for E° below 0.5 MeV was quite different from the 
experimental results. The also obtained the variation of 
the alpha energy distribution with R and compared them with 
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the results of Fraenkel (1967). The effect of assuming that 
the alpha particle is created off the fission axis ( ^3f ) s 
on these results was shown to he very small. It was also 
not iced that the value of E° was influenced by the choice of 
X , and based on prompt neutron data, the value of 0.5 MeV 
was found to be appropriate and the point of emission (for 
R = 1.4) was closer to the heavy fragment (X = 12f measured 
from the centre of the heavy fragment; the zero field point 
would be cl t about 15f ) . 

Cavallari et al. , (1971) used a range of values for 
the initial parameters, E°, X , Y q and Q°. These ranges 
were divided into a number of intervals. The trajectories 
were computed for fixed R using different sets of initial 
values selected from these intervals and for 5 different 
values of Ey. Using experimental results (specifically E ) 
they assigned a ’weight' to each trajectory and were thus 
able to obtain statistical distributions for the initial 
parameters. These distributions could be obtained only when 
3° = 0 MeV. The B° distribution showed a strong peak at 
0. 5 MeV and a weak one around 2 MeV. The X Q distribution 
indicated that the alpha particle was emitted at the zero 
field point. The ©^ distribution was double peaked. The 
most probable values obtained by them are given in Table 5.1. 

More recently Choudhury (1976) has been able to 
reproduce the experimental distributions as well as some 



Table 5.1 

Initial values of calculations supporting the Statistical model 
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observed correlations using zero initial kinetic energy for 
both the fragments and alpha particle. In particular, the 
observed anticorrelation between the average final fragment 
kinetic energy and the final alpha energy has also been 
reproduced. 

The results of the calculations which lend support 
to the predictions of the statistical theory are summarized 
in Table 5.1. 

5*5*2 Calculations supp o rting Adiabatic models 

Halpern (1963) developed the three-point-charge 
model and computed the energy and angular distribution of 
the alpha particle in ternary fission of He est ima ted 

the initial alpha energy to be 4.4 MeV from the volume of 
the neck at scission using the uncertainty principle. The 
initial angular distribution of the alpha particle was 
assumed to be isotropic and the final angular distribution 
was found to be insensitive to the location of the point of 
emission of the alpha particle and to its initial energy. 

The average final alpha energy showed appreciable variation 
with initial alpha and fragment energies. In order to 
reproduce the observed distribution in alpha energy he 
concluded that the initial fragment energy had to be about 
15 to 20 % of the final fragment energy if the initial alpha 
energy was 4.4 MeV. 
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Boneh et al., (1967) computed trajectories of the 

252 

alpha particle in ternary fission of Cf using the three- 
point-charge model in the two dimensional approximation. A 
'standard' set of values were chosen for the seven initial 
parameters. This 'standard 1 set xvas obtained from some 
preliminary calculations. The value of R was chosen to be 
1.4 since this is roughly the mean value observed experimentally. 
The preliminary calculations consisted of the determination 
of the effect of the variation of one initial parameter on 
the final values, the other initial values being fixed. 

Using this standard set they reproduced most of the experi- 
mental results of Eraenkel (1967). In order to reproduce 
the correct widths in angular distribution and the anti- 
correlation between the final fragment energy and the final 
alpha energy, they concluded that the initial energies had 
to be high, that is, E° 3 MeV and Ep 35 MeV. The 
emission angle of the alpha particle was found to depend 
principally on the emission point (on the fission axis) and 
the mass ratio. They had also attempted to calculate the 
initial distributions using empirical relations connecting 
various initial and final parameters. Using some of these 
distributions they could successfully reproduce the experi- 
mental angular and energy distributions of the alpha particle. 

The initial fragment energy distribution was found to be 
Gaussian with a mean of ^39 MeY and the initial alpha energy 
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distribution was a Maxwellian dependent on the point of 

emission of the alpha particle. Mo unique distribution for 

X could be obtained, 
o 

Katase (1963) carried out similar calculations for 

the fission of and ^"Cf. The results of the calcula- 

2 3 G 

tions for U were compared with experimental results and 

the dependence of certain final parameters on initial ones 

were identified. On this basis empirical relationships 

between the final and initial angles and between the alpha 

energy and various initial parameters were obtained. The 

parameters of the latter expression (about 8 of them) were 

obtained by a least squares method, that is, the values of 

the coefficients that gave the best fit to the values of 

obtained from trajectory calculations, were determined. 

252 

Using these expressions the final distributions for Cf 
were determined and compared with experiment. The agreement 
with experiment was good and he concluded that the fragments 
had about one third of their final energy at scission and 
the initial alpha energy was about 3 MeV. 

Raisbeck and Thomas (1968) obtained a good fit to 

the experimental kinetic energy spectrum of the alpha 

252 

particles from Cf in terms of a model in which the alpha 
particle materialises at a point (the zero field point) on 
the fission axis, at a time * t* after binary fission has 
occured. At the instant of scission the fragments were 



101 


assumed to have zero velocity and the alpha energy was 

described by a Gaussian distribution. They obtained the 

best fit with t = 0.4xl0 -2 "'" sec and E°~'2 MeV. At the instant 

a 

of emission of the alpha particle the fragments had a mean 
kinetic energy of ~3 MeV. They had also carried out similar 
calculations for other light particles such as protons, 
tritons and °He in an attempt to determine whether a single 
set of parameters could account for the energy spectra of 
these pa. rticles. -They concluded that the emission mechanism 
of all light particles was essentially identical. 

Block! and Krogulski (1968), carried out computa- 

236 

tions for alpha emission in the fission of U. Using the 
sudden approximation i.e. abrupt scission they assumed the 
fission fragments to be at rest at scission. The initial 
alpha energy was described by a Maxwellian distribution of 
the type PJ exp (-E q /T) with T = 0.35 MeV. They could 
reproduce the experimental distributions quite well. The 
calculations were also used to study the emission of other 
light particles such as 2 H, ^H, He and other heavier 

particles upto "^0. They concluded that the emission 
conditions - of these nuclei were very similar to those of 
the alpha particle. In a subsequent calculation (Krogulski 
and Blocki, 1970) for 252 Cf, the effects of the various 
initial conditions on the final energy and angle were 

The results led them to conclude that the initial 


examined 
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conditions for the emission of alpha particle, tritons and 
He were quite similar. However, an extremely elongated 
configuration was needed for the emission of protons and 
deut rons ; which led them to conclude that the emission 
mechanism for these particle was different from that of alphas 
for eg. emission from the surface of the fragments instead 
of from the neck. 

G-azit et al. , (1971) used trajectory calculations 

236 

to complement their measurements on U fission. Their 
calculations included corrections for alpha recoil effects. 
They obtained good agreement with their experimental results. 

Musgrove (1971) obtained a set of initial values 
quite close to those of Boneh et al. , The angular distri- 
butions he obtained were however, narrower. 

Rajagopalan and Thomas (1972) using initial 
conditions similar to those used by Raisbeck and Thomas 
(1968) obtained angular distributions which agreed very well 
with their measurements on angular distribution. 

The results of the calculations which lend support 
to the dynamical theories are given in Table 5.2. 



Initial values of calculations supporting Adiabatic models. 
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5 • 6 Present calculations 
5.6.1 Motivation 

The common feature of most of the calculations that 
were reviewed in the previous section, is that of assigning 
a set of fixed values to the initial parameters. Using this 
fixed set, certain final distributions and correlations have 
been obtained. In some of these calculations eg. Halpern 
(1963) and Vitta (1970) the initial set was chosen based on 
the assumptions or predictions of some model of fission. 

The initial set of values used in calculations 

supporting the predictions of either the statistical or 
adiabatic models are quite different, yet they both repro 
duce final distributions and some correlations quite 
satisfactorily. The inability to resolve the apparent 
ambiguity in the initial values could be due to the imposed 
limitation of choosing essentially fixed values for the 
parameters, whereas in reality each parameter would have a 
distribution, however small. This argument is also supported 
by the fact that in many calculations it was found necessary 
to assume some small widths for certain initial parametero 
in order to obtain the correct widths for the final distri 
butions eg. Vitta (1970) had to assume a spread in of 
about 2f to account for the observed spread in similarly 

Boneh et al. , (1967) had to assume certain widths for E a 
account for the width in © T . The mean values of these 
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distributions - usually Maxwellian or Gaussian - represent 
the fixed values. While reproducing the experimentally 
observed anticorrelation between Ep and E , Boneh et al„, 
and others used fixed values for the paiameters; however, 
they had remarked that distributions in initial parameters 
would effect their results. Thus such methods of obtaining 
final correlations are open to question. 

It would therefore appear realistic to begin with 
a distribution of values for the initial parameters and 
then attempt to reproduce the final distributions and 
correlations. The nature of initial distributions obtained 
this way can then possibly indicate the validity of one of 
the two models discussed above. 

So the method of attack should be to remove all 
artificial restrictions on the values of the initial 
parameters and make a model independent calculation. Such 
a procedure had been adopted by Cavallari et al. , (1971). 

The philosophy of the present calculations is similar. T^.e 
initial parameters are allowed a range of values covering 
ail regions studied so far. Each value in the range initially 
has an equal probability of occurance. furthermore the 
parameters are initially uncorrelated. These assumptions , 
in fact represent one extreme choice of initial values, 
however, they impose no restrictions on the behaviour of 
the calculated final values. 
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5.6.2 Choic e of values 

The range of values for each of the seven parameters, 

were essentially decided by performing some preliminary 

calculations (Gadgil, 1973) similar to those done by 

Boneh et al. , (1967). In these calculations all initial 

parameters were held fixed except one and the sensitivity of 

the asymptotic values of E , E , and G T etc. to this para- 

o: p n 

meter were checked. A representative set of results of • 
these calculations is shown in Pigs. 5.2 (a) to (f). The 
(fixed) values of all parameters except the one being varied 
are also indicated. It was seen [Figs. 5.2 (e) and (f)]that 
for variations in Q? between 40° and 140° and Y between 

ij o 

0 and 2.5 f; the final values in E^, Ep and Gp did not change 
as significantly as for the variations in other initial 
parameters. Hence in the subsequent calculations Op was 
allowed an isotropic variation between 45° and 135° and Y 

o 

was fixed at 0 f i.e. the alpha particle was assumed to be 
emitted from the fission axis. The assumption on Y q 
contradicts the philosophy of allowing all parameters free 
variation within a range of values, however, since the 
preliminary calculations and many previous calculations 
indicated that for Y q = 0 to 5 f the results were not 
c?„itered significantly and because of limitations imposed on 
computing time; such an assumption had to be accepted. 
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The other parameters, namely D, E°, E^ and X Q were 
allowed values in the ranges given in Table 5.3. Each range 
was divided into ten intervals thus allowing each parameter 
to take any of the ten values. The calculations were perfo- 
rmed for seven different values of mass ratio (Jlrishnara julu 
and Mehta, 1975). 


Table 5.3 



The input parameters 

and their 

range of values 


Parameter Symbol 

Range of values 

1. 

Mass ratio 

R- 

1.0 to 1.8 

2. 

Interfragment distance 

D 

17f to 28 f 

3. 

Initial distance of alpha 
particle from zero field 

X 

-7. Of to 6.5 f 


point 

0 


4. 

Initial distance of alpha 

Y 

0 f (fixed) 


particle from fission axis 

0 


5. 

Initial alpha energy 

E° 

a 

0.4 to 4.0 MeV 

6. 

Initial fragment energy 

E° 

"F 

0.5 to 40 MeV 

7. 

Initial angle of emission 

©° 

W T, 

35° to 145° 


of alpha particle 

JJ 


5.6 

.3 Details of calculations 




Each trajectory was computed using a set of initial 
values for the parameters, characterized by a five digit 


\ 
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number - corresponding to the five varying initial parameters 
R, E°, Ep, D and X Q (©^ varied isotropically) . It took 
0.44 sec to calculate a trajectory on an IBM7044 computer and 
70,000 such trajectories were calculated. The five digit 
number specifying the initial values of the trajectory and 
the final values of E^, and 0^ were written on magnetic 
tape and stored thus enabling subsequent analysis to be 
performed with relative ease and minimal computer time. 

In order to pick out the physically significant 
trajectories, the final results of each trajectory were 
compared with experimental distributions and a relative 
•weight' was assigned to the trajectory as follows. The 
experimental distributions in and Ej, for different values 
of mass ratio (Mehta et al. , 1973) were used. These 
distributions were assumed to be Gaussian and were defined 
by specifying the mean value and the variance of the distri- 
butions for various mass ratios. The ‘weight' was a product 
of terms of the type 

= exp ( — | x— x [ /2a ) 

- 2 

where x and a are the mean value and variance of the 
experimental distribution being compared with the calcu- 
lated value x of the variable (E^ or Ep in the present case). 
All trajectories with P^O.l where P = P P , were used in 

x j 

subsequent analysis. 
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Results 


5.7.1 E and Ep distributions 

In order to check the calculations the mean values 

of the E and E-r, distributions obtained from the calculations 
a E 

were compared with the experimental ones used in assigning 
weights. The comparison is shown in Table 5.4. 

Table 5.4 


Comparison of experimental and calculated values of E^, Ep 


R 

E 

a 




Calc 


Calc 

1.0 

17.5+0.5 

17.54+1.03 

173 +1.0 

172.85+1.32 

1.2 

16.0+0.1 

16.2+0.98 

173.5+0.5 

173.29+1.26 

1.3 

; 15.95+0.1 

16.12+0.98 ! 

170.5+0.5 

170.39+1.24 

1.4 

15.95+0.15 

16.16+0.98 

166.75+0.5 

166.87+1.18 

1.5 

16.0+0.2 

16.43+0.98 

162.75+C. 5 

163.99+1.08 

1.6 

16.05+0.25 

16.98+0.88 

158.75+1.0 

161.70+0.94 

1.8* 

16.25+0.55 

16.94+0.73 

152.5+1.0 

162.29+1.15 


Statistics were poor for this mass ratio. 

As can be seen from the table the agreement is quite 
satisfactory. 
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5.7.2 Variation of ©^ with R and the dependence of E a on ©^ 

The angular distributions of alpha particles for 
different values of mass ratio were obtained and were compared 
with the experimental results of i’luss et al. , (1973). These 
results, for a few values of R, are shown in Fig. 5.3. The 
agreement is satisfactory indicating the validity of the 
method used. The variation of the average ang3.e of emission, 
©£, with mass ratio, R, is shown in Fig. 5.4. The results 
are compared with those of Fraenkel (1967) and Fluss et al. , 
(1973). The average angle of emission decreases with 
increasing mass ratio. This trend of variation is in agree- 
ment with experimental results and for R = 1.3 amd 1.4 the 
actual values are very nearly equal to the observed ones. 
However, there is disagreement between the calculated and 
experimental results in the region 1.5<R<'2 and also for 
R = 1.0. This could be due to the inherent experimental 
inaccuracies in these regions of mass ratio. However, the 
value of ©^ 90° obtained for R = 1 in the present 

calculations represents a physically meaningful result since 
both fragments exert an equal force on the alpha particle it 
should come out in a direction perperndicular to the fission 
axis. If one takes into account the recoil due to the 
emission of the alpha particle - the departure from colinearity 
is about 4.3° (Fluss et al. , 1973) - then the value of ©^ 
for R = 1 could be about 92°. In the case of large values of 
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I he variation of tho average angle of emission. of the alph 
particle (Q^) with, mass ratio (R) . Also shown are ti, 

results of Fraenkel (1967) and Fluss et al., (1973). 
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R, Fraenkel's results (1967) show an appreciable decrease 
in whereas those of Fluss et al., (1973) indicate a very 
small decrease. This is because of the stated poor mass 
resolution in the latter case. The present results lie 
between these two experimental results. 

A second correlation which is significant is that 
between the angle of emission of the alpha particle and its 
energy. The variation of the average alpha energy, E , 
with the angle of emission is shown in Pig. 5.5. There 
is quantitative agreement between the present calculated 
results and the experimental ones. The results of the 
calculations of Boneh et al. , (1967) are also shown for 
comparison. It is seen that the average value of the alpha 
energy is minimum at the most probable angle of emission 
(~'84°) and increases as changes in either direction. The 
calculated values of are, in general, higher than the 
experimental ones, especially when compared with the results 
of Fluss et al. , (1973). This could be because the alpha 
energy averaged over all angles obtained by Fluss et al. , 
is about 1 MeV less than the accepted value of 16 MeV. 

It has been thus demonstrated that all these 
experimental results can be reproduced with the initial 
parameters having distributions, using the method adopted 
herej that is, of assigning statistical 'weights ’ to the 
trajectories using a few experimental distributions namely 
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the E and distributions . When only the E and © T 
distributions were used it was seen that the average total 
fragment energy E^, was independent of mass ratio contrary 
to experimental results. This result was interpreted to 
mean that in order to account for the deformation energy 
of the fragments (which is inherentljr neglected due to the 
point 'charge approximation) it was necessary to include the 
Ep distributions while assigning weights. 

5.7.3 The anticorrelation between Ej, and E 

The total kinetic energy of the three-particle system 
is a conserved quantity in the framework of such trajectory 
calculations. If the average initial -energy of the fragments 
is uncorrelated with the initial kinetic energy of the alpha 
particle then an increase in one of these quantities during 
motion must be accompained by a decrease in the other 
quantity. If the initial energies are small and fixed 
(that is, the scission configuration were always the same) 
then an increase in alpha energy must be compensated by a 
decrease of equal size in fragment energy i.e. 6Ep/6Ea=-l. 
However, if the initial energies were large or if they had 
a distribution of values (implying different scission 
configurations) then the degree of correlation would, in 
general, be smaller than unity.' 
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Experimental measurements on this anticorrelation 

(Fraenkel, 1967 j Mehta et al. , 1973) have shown that its 

252 

value averaged over all mass, ratios for Cf is -0.44. Boneh 
et al. , (1967) had in their calculations, obtained the correct 
anticorrelation by keeping all initial variables fixed except 
E° which was varied between 1 and 5 MeV. With lower values 
for the initial energies they obtained a stronger anti- 
correlation. This led them to conclude that at the moment 
of scission the three particles had already acquired about 
20/£ of their final kinetic energy. Choudhury (1976) has 
obtained fair agreement with experimental results starting 
with zero kinetic energies and narrow distributions for 
D, X and Y . He attributed the stronger anticorrelation 
obtained by Boneh et al. , (1967) to their fixing of D and Y . 

It must be emphasized that in the results quoted above the 
initial fragment and alpha energies were assumed to be 
uncorrelated. The presence of distributions in the initial 
parameters and/ or correlations between them could alter the 
conclusions regarding the scission configuration. Therefore, 
in order to get a clearer picture of the scission configura- 
tion it would be necessary to examine the effect, of distri- 
butions in the initial parameters and also, if possible, the 
correlations between them. 

The method adopted in. the present calculations has 
therefore been directed by such requirements. The initial 
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parameters have been allowed a distribution in values and 
the anticorrelation obtained is shown in curve A of Pig. 5.6. 
The degree of anticorrelation obtained is very small. 

However, the correct anticorrelation was obtained when the 
fragment trepanation D, was held fixed (curve B) while all 
the other parameters were allowed to have distributed 
values. In particular, no regard was given to the initial 
energies being high or low. This result indicates that 
when the interfragment distance D is held fixed a negative 
correlation between and E is obtained with no restric- 

.r QC 

tions on the E° or Ep distributions. One, therefore, 
cannot conclude that at scission the particles have very 
small energies or that their energies are about 20j^of 
their final values. 

The result that the anticorrelation almost disappears 
when all the parameters, including D, are allowed to have 
distributed values; could be due to an initial positive 
correlation between E° and Ep. Boneh et al., (1967) had 
pointed out that a distribution in the values of D could 
give rise to a positive correlation between E° and Ep. They 
had, however, stated that this positive correlation would 
be partially compensated by the negative correlation due to 
the distribution in initial heavy fragment velocity (i.e. 
in effect a distribution in E°) . There exists no quantitative 
support for this conjecture. The present calculations 




The antxcor relation be tween the average total fragment energy {£ 
and the alpha amerg y j for a fixed value of D and for D having 
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indicate that the distributions in the other parameters 
(apart from D) could account for this (possible) positive 
correlation. 

In the results presented above only the effects of 
distributions have been considered however one cannot exclude 
the possibility of some initial correlations which could be 
responsible for the final anticorrelation. With fixed D it 
was possible to obtain the correct anticorrelation hence this 
could mean that fixing D gives rise to a negative correlation 
between and E^ which more than compensates the positive 
correlation in these two quantities due to the distribution 
in all other parameters. In the foregoing conjecture it has 
been assumed that the initial correlation between E° and E° 
largely determines the sign of the final correlation between 
Ep and E^. 

As illustrated in Eig. 5.6 the correct anticorrela- 
tion was obtained with D = 25.8+0.2 f. When the variation 
was increased to 1.0 f the anticorrelation became very small 
indicating that the anticorrelation depends quite sensitively 
on the spread in value of D. In other words, the D distribu- 
tion necessary to produce the correct anticorrelation is 
very narrow. This result indicates that the alpha particle 
is emitted when the fragments effectively maintain a 
constant separation. This could be pictured as follows t 
the motion of the fragments just after scission and the 
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change in their deformations in this period, could result 
in zero nett displacement of the charge centres of the 
fragments thereby resulting in a fixed D during this small 
time interval . 

The variation of the value of the anticorrelation with 
mass ratio was studied. The value of D chosen for this 
purpose was 26.2+0.2 f. The results, shown in Pig. 5.7 are 
compared with the experimental results of Mehta et al. , (1973) 
The calculated values agree quite well with that of experiment 

5.7.4 Initial distributions 

Bach trajectory was assigned a 'weight' by comparing 

the final B and E-™ values with experimental distributions, 
a £ 

This 'weight' , denoted by P, with P = P x ?y as explained 
earlier, expressed the probability that the initial values 
characterizing the trajectory had of reproducing the mean 
values of the experimental distributions in B^ and Bp. In 
choosing the trajectories for analysis a lower limit on P 
was imposed; only those trajectories which had P>0.1 were 
considered. This cut off though not physically significant 
as such, does however, imply that only those initial values 
having more than 30^ chance of resulting in each of the mean 
experimental values of E^ and Ep were considered. This cut 
off was also chosen with a view to obtain good statistics. 

The mean values of and Ep were not. altered when P 0.01 
or P> 0.001. 
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These initial values were then plotted against 
their frequency of occurrence, that is LT(x) vs x, to give 
the initial distributions. Thus, while counting the number 
of times a particular initial value o c cured , no regard was 
given to P. That is, once the trajectories had been 
selected they were treated on equal footing. In this manner 
the distributions in E°, E^, X Q and D were obtained. 

The distribution in E° is shown in curve A of 
Pig. 5.8 (a) (curve B will be discussed later). It is clear 
that this distribution does not give the much-sought— for 
information on whether the initial energy is low . or is an 
appreciable fraction of the final value. The indication is 
that the distribution is broad with possible peaking in the 
1 and 2 MeV regions. The E° distribution, shown in curve A 
of Pig. 5.8 (b) is also quite broad and shows possible 
multiple peaking. In order to examine this possibility the 
variation in E^ was extended to 60 MeV and the result is 
shown by the dotted extension in the figure. 

The distribution in the point of emission of the 
alpha particle along the fission axis is shown in Fig. 5.8(c) 
It is seen that the most probable initial starting position 
for the alpha particle is close to the zero-field point. 
Several previous trajectory calculations eg. Boneh et al., 
(1967), Raisbeck and Thomas (1968) also gave the same result. 
Pong (1970) remarked that the main difference, which accounts 
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for different values for E° between other caicuoations and 

a 

his, is the choice of X . The zero field point, as the 
most probable emission point, corresponds to a position 
closer to the light fragment j Pong however, requires it to 
be nearer the heavy fragment. 

The distribution in initial values of interfragment 
distance, D, shown in Pig. 5.8 (d) is broad with a mean 
value of '-'24 f and a width of about 6 f. This value is 
quite close to the one obtained by Pong and others. Boneh 
et al. , obtained 26 f while Raisbeck and Thomas obtained 21f . 

These distributions were first checked to see whether 
they reproduced the final distributions in E a and E-pj the 
ones used in assigning 'weights’ to the trajectories (see 
Table 5.4). Then they were used to predict angular 
distributions and a few final correlations. The calculated 
distributions agreed quite well with experimental distribu- 
tions as elaborated in section 5.7.2. However, the observed 
anticorrelation between E™ and E could not be reproduced 

£ (X 

by these distributions as explained in the last section. 

In order to obtain the initial distributions which 
reproduce the observed anticorrelation it was thought fit to 
use the anticorrelation condition as a criterion in assigning 
'weights’ to the trajectories in addition to the E a and Ep 
distributions. The values of the anticorrelation for 
different values of R were taken from the experimental 
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results of Mehta et al. , (1973) and were used in the same 

manner as was done with the E and E-r, distributions. 

a jj 

The initial distributions obtained when the anti- 
correlation condition was also used, are shown by curves B 
in Pigs. 5.8 (a) to (d). Here again there is an indication 
of two peaks both in the E° and distributions with some 
changes in their positions and relative amplitudes. By 
and large there is no significant change in these distribu- 
tions. The X Q distribution too does not change when the 
anticorrelation condition is used in assigning weights. 

This implies that the anticorrelation condition has no 
significant effect on these initial distributions. 

The distribution in initial values of D when the 
anticorrelation condition is used is shown in curve B of 
Pig. 5.8 (d). The distribution is broad, however, there 
are indications of two peaks one around 23.5 f and the other 
around 25.5 f. It was seen that with fixed D it was possible 
to get the correct anticorrelation, however, using the 
anticorrelation condition the distribution in D is not 
narrow and single peaked as one might expect. The distri- 
bution in D for different values of R indicates that the 
(broad) peak shifts (~1 to 2 f) towards lower D values 
for mass ratios different from R = 1.4. Such a result could 
possibly imply that in symmetric fission the deformation 
(elongation of the neck) is relatively not large and the 
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increase in "the total number of prompt neutrons for symmetric 
mass division (Terrel curve) is probably due to increased 
fragment excitation coming from viscous effects. 

In order to throw more light on this apparent paradox 
(ei.e. D not being single peaked) the initial distributions, 
with fixed D, were studied. The values of D chosen for this 
purpose were 23.25+.25f and 25.5+. 25 f, corresponding to the 
approximate peak positions referred to in the last paragraph. 
The distributions in E° and Z Q are shown in Figs. 5.9(a) and 
(b) . The Ep distributions for these two cases were narrow 
and had mean values of 15 MeV for the lower D value and 
35 MeV for the higher D value with corresponding mean alpha 
energies of 1.0 and 3 MeV respectively. The Z Q distributions 
in both cases are broad with mean values corresponding to a 
position closer to the light fragment in the case of the 
lower D value and to a position closer to the heavy fragment 
for the higher D value. In the latter case, however, there 
is an indication of two peaks on either side of the zero- field 
point with the peak nearer the heavy fragment being more 
prominent. The distribution reflects this double 
peaking (Pig. 5.9(c)). 

The significant aspect of the results shown in 
Figs. 5.9(a) and (b) is that they indicate that both classes 
of initial values i.e. those of Boneh et al. , (1967) 

(adiabatic model) and of Fong (1970) (statistical model) 
can be obtained. Boneh et al. , used initial alpha and 
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fragment energies of 3 Me? and 35 Me? with D^26f and Fong 
had used initial energies ^0.5 and 1 Me? with D — 24f. 

Boneh et al. , had obtained the zero field point as the most 
probable emission point while Fong had required this point 
to be closer to the heav 3 r fragment. This latter result does 
not agree with the most probable X Q obtained in the present 
calculations for Ds^23.5f. 

Gavallari et al. , (1970) had also obtained distribu- 
tions in the initial parameters and their average values (see 
section 5.5.1). They had however used three average final 
alpha energies, of 9, 15 and 18 Me? in assigning weights to 
the trajectories. The other final values eg. for Ej, and ©^ 
were just used to select those trajectories which gave the 
best approximation to the experimental values (Ej, = 177.3 Me? 

©^ = 92+30°). Hence* at best the initial distribution obtained 
are those which can reproduce any of the three final values 
of alpha energy and approximately the final fragment energy 
and final angle. Using the initial values they obtained 
(see Table 5.1) it would not be possible to obtain the 
correct anticorrelation * in fact the value of tne anticorre- 
lation obtained for very small initial energies would bess -1 as 
explained in section 5.7.3. They did not check their 
distributions against the final fragment distributions 
(except for the requirement that 1 . 32 ) or the 

anticorrelation condition as has been done in the present 


calculations. 
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To summarize the results on initial distributions 
one firstly notices that the initial distributions are 
broad , covering both the predictions of the adiabatic and 
statistical model as elaborated above. The noteworthy fact 
about them is that even when the anticorrelation condition 
was imposed as a criterion, these distributions remained 
unaltered. This fact suggests that the initial distributions 
are in fact broad and choosing either low or high initial 
energies within the distributions would not alter the final 
values substantially. It would therefore be incorrect to 
say either that the particles have attained an appreciable 
fraction of their final energy or that they are almost at 
rest at scission. 

Most calculations including the present are support 
the result that the most probable emission point for the 
alpha particle is the zero-field point. However, there is 
some disagreement in the value of the interfragment distance 
at scission. The present calculations indicate that in 
order to reproduce the final anti correlation a very narrow 
distribution of interfragment distance is required and its 
mean value is close to that obtained by some calculations 
eg. those of Boneh et al. , (1967). 

5.7.5 Correlations between initial parameters 

One of the simplifying assumptions made in all 
trajectory calculations including the present one, is that 
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of treating the initial parameters as uncorrelated quantities 
within the constraints of the conservation laws. This assump- 
tion is made because of the lack of apriori knowledge about 
such correlations. One thus begins with uncorrelated 
initial parameters. However, in order to explain some of 
the final correlations, some ' qualitative conjectures about 
initial correlations have often been made (eg. Boneh et al. , 
1967). As yet, none of the calculations have succeeded in 
coming up with some definite information on initial 
correlations . 

In the last section, the effect of distributions in 
initial parameters on the anticorrelation was examined. In 
order to got more information on the final anticorrelation 
it was thought worthwhile to examine the effect of some 
possible correlations in initial parameters. Intuitively, 
one would imagine they an initial negative correlation 
between E° and E-2 would lead to the final anticorrelation. 

Ci r 

However, correlations between the other parameters could 
also influence the final anticorrelation. 

In order to study the (possible) initial correlation 
between E° and Ep the distribution in Ep when S° 1 MeV and 
for > 1 MeV were obtained. These distributions, shown 
in Eig. 5.10, indicate a positive correlation between E° 
and Ep. The distributions remain unchanged when the anti- 
correlation condition is also imposed. The mean values of 
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Fig. 5.10. The correlation between the initial alpha 

energy (Eg) and the initial fragment energy 
(gp) for R = 1.4. The ordinate is in arbitrary 
units. Case (I) and case (II) as in Fig. 5.8, 
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from these distributions are close to the ones obtained by 
Bong (1970) and others (E°^l MeV) and by Boneh et al. , ( 1967 ) 
and Raisbeck and Thomas (1968) (E° > 1 MeV). This result 
indicates that although there is an initial positive correla- 
tion between E and Ep the final values are negatively 
correlated. This initial positive correlation is probably 
compensated for by some other correlation which is strongly 
negative. In terms' of the energy transfer from deformation 
to kinetic energy, this result indicates . that the extent to 
which this transfer takes place is either small, resulting in 
small E° and Ep ; or is large resulting in large E° and 
Ep. Both situations seem equally valid and lead to the 
final anticorrelation. 

As mentioned in the previous section the initial 
positive correlation between E° and Ep remains essentially 
unchanged when D is held fixed. Eor the two values of D 
studied, namely 23. 5f and 25.5f» the average values of E° 
and Ep were 1.0 and 15 MeV and 3.0 and 35 MeV respectively. 
These results may be pictured in the following manner; when 
the transfer from compound nucleus excitation / deformation 
to kinetic energy is small the kinetic energy of the nascent 
fragments is small with a relatively small interfragment 
separation and when the energy transfer is large implying 
a larger kinetic energy of the nascent fragments the separa- 
tion is relatively larger. The results of Bong (1970) and 
others indicated a lower value of D than that obtained by 
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Boneh et al. , (1967). 

The X distribution for these two regions of E° is 
0 ^ 

peaked around the zero field point to within 0.5f. The X Q 
distribution is therefore not affected by a choice of 
either low or high initial energies contrary to the sugges- 
tion by Eoneh et al., (1967) that the initial alpha energy 
is dependent on X Q . These distributions did not change when 
the anticorrelation condition was imposed indicating their 
valid ity . 

Another possible initial correlation is that between 
X Q and 1. It was seen that in order to obtain the correct 
anticorrelation D had to be held fixed and the X Q distribu- 
tion was broad with possibly two peaks (fig. 5.9(b)). When 
the anticorrelation condition was used as a criterion in 
assigning weights, the X Q distribution was peaked around the 
zero-field point whereas the h distribution indicated the 
possibility of two peaks (fig. 5.8(d)). A similar result 

i.e. two values for X for a fixed D, has been obtained by 

o 

Choudhury (1976). Some empirical correlations between X Q 
and D were tried but they failed to produce the correct 
anticorrelation. The exact nature of this correlation is 
not clear but it does seem that in order to produce final 
energies satisfying the anticorrelation condition the 
potential energy configuration at the moment of alpha emission 
(determined by D and X Q ) is important. 
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5 . 8 The e ffect of fragment deforma t ion on initia l dis tributions 

In section 5.2 certain approximations used in trajectory 
calculations were explained. These approximations were 
necessary in order to simplify the problems in calculation by 
reducing the number of parameters. If the fragments are 
considered to be undeformed spheres, then their charge centres 
would coincide with their centres and hence could be considered 
as point charges. This means that higher order multipoles 
arising out of a non-spherical charge distribution are 
neglected. The non- inclusion of dipole and quadrupole inter- 
actions would reduce the potential by 1% and 4 ^respectively 
(see section 5.2). Since 3 and E- are about 16 MeV and 

OC r 

180 MeV respectively, the errors due to this approximation 
could amount to about 1 MeV in E„ and about 8 MeV in E-. 

These figures are actually overestimates since after scission 
both D and the deformation of the fragments are rapidly 
changing and the potentials due to these multipoles would 
become insignificant after a small time interval i.e. a few 
iterations in the calculations. It was seen (Table 5.4) that 
the average final energies obtained from these calculations 
agreed very well with those of experiment. Hence, as was 
also pointed out by Halpern (1963) , the mean final energies 
would not be significantly altered by the inclusion or non- 
inclusion of higher order multipoles in the interaction. 
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However, the inclusion of deformations could effect 

the distribution in D and since the charge centres and 

o 

consequently the zero field point would shift depending upon 
the extent of deformation. Hence it would be instructive to 
compare the initial distributions obtained using deformed 
fragments with the earlier ones. 

Geilikman and Khlebnikov (1964) had carried out 
trajectory calculations using deformed fragments. They had 
considered quadrupole and octupole deformations and had 
neglected the change of these deformations with time. Their 
results on energy and angular distributions compare very 
well with those of Halpern (1963) who used the point-charge 
approximation. 

The aim of the present calculation was to determine the 
effect the inclusion of deformations had on the initial 
distributions. For a mass ratio R = 1,4 it was assumed that 
only the light fragment was deformed into a spheriod with 
e = 0.7 and that the heavy fragment was spherical. This 
case is only representative and is meant to determine whether 
the inclusion of deformation had any effect at all on the 
initial disteibutions . It was further assumed that the 
deformation was independent of time i.e. the fragments 
retained their shapes. 

The potential due to a spheroidal nucleus at (r,Q) is 
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where Ze is the charge of the nucleus and 'a' and 'b' denote 
the semi-major and semi-minor axis respectively. By volume 
conservation we have 
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E(r, S ) = -3*2 [i^i 
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Since P 2 (Cos ©) = 1/2 (3Cos 2 ©-l); P£(Cos ©) = 3 Sin © Cos © 


P 4 (Cos ©) = 1/8 (35Cos 4 «-30Cos 2 ©+3) ; 


P^ (Cos 0) = -5/2 Sine © (7Cos 3 ©-3Cos ©) 


md r = Sin © Cos fb 1 + Sin © Sin /) j + Cos © k 


© = Cos © Cos /> i + Cos © Sin jb j - Sin © k 


where i , j and k are unit vectors along the x,y and z directions. 
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56(r/R Q ) 4 r 4 r 2 
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3Ze r z , e 2 
^ ^ 10(r/R o )^ 



(x 2 +y 2 )-z)+ j 

56(r/R Q ) 4 



30z' 


+ 3z 


- (is! . ^ })] 


r -r' 

These relations were used in the calculation of forces on the 
three particles. 


The values of the initial parameters used were the 
same as those used earlier. Ten thousand trajectories were 
computed for R = 1.4 and each trajectory was assigned a 
'weight' as before. The maximum percentage error in energy 
conservation was ~ 2%, 

The initial distributions obtained are shown in Rig. 
5.11(a) to (d). One can at once see that there is no 
significant change in the distributions as a result of 
including deformations. These distributions also remain 
unaltered when the anticorrelation condition is imposed. The 
initial positive correlation was also reproduced. The final 
Q distribution did not show any significant change though the 
mean value was slightly reduced from 86° to 84.6°. 

The method adopted here to examine the effects of 
deformation of fragments has not produced any new result 
except to indicate perhaps that unless the effects of defor- 
mation are larger there would be no significant change in 
the initial distributions. The magnitude of the interaction 
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between the fragments has increased due to the inclusion of 
higher multipoles, however this increase docs not alter the 
the results. In order to study the effects of deformations 
more thoroughly, time dependent deformations of both fragments 


should be included 



CHAPTER VI 


SUMMARY OP RESULTS 


Channel effects have been observed in prompt neutron 

emission, average kinetic energy of fission fragments etc., 

235 

in slow neutron (thermal to 1 MeV) induced fission of U. 

The aim of the Present work was to determine whether such 

effects could be observed in the yield and average energy of 

235 

light charged particles (LCP) emitted in the fission of U 
23Q 

and ' vfu in this neutron energy region. 

The LCP yield shows definite structure outside the 
limits of error for both the nuclei studied. The LRA yield 
shows distinct structure around 200 keV neutron energy — an 
increase of about 20 / 0 above the thermal yield in the case of 
2,5 U and about 1 % in the case of 259 Pu. The 259 Pu results 
also indicate increases around 600 keV and ,1 MeV and a drop 
("vl4%) between 250 and 300 keV. The variation of average 
LRA energy with neutron energy shows some structure (~ 0.6 MeV) 
around 200 keV. With better neutron energy resolution 
it would be possible to examine this structure more thoroughly 
and draw definite conclusions about variation of average 
energy. The width, on the other hand, seems to be insensitive 
to variation in neutron energy. 

235 

The increase in LRA yield around 200 keV for U 
can be understood in terms of the increase in fission via 



148 


even parity states preferentially populated due to the p— wave 

interaction in this energy region. Since the positive parity 

states lie lower than the negative parity states the 'extra* 

energy made available goes to enhance the ISA emission 

probability. The drop in LBA yield between 250 and 300 keV 

239 

in the case of ' Pu can also be explained in terms of the 
relative predominance of the p-wave interaction. In this 
case fission takes place predominantly via negative parity 
states (eg. -1”) which lie higher than the positive parity 
ground state (0 + ), thus the 'decrement' in energy ^is refle- 
cted in the reduction in the LEA emission probability. 

The extracted 'triton' yields also indicate 

structure in the neutron energy region investigated. The 

235 

strucutre is more pronounced in U fission as compared to 

25 %tL fission. The increase ( 2 to 3 times) around 500 keV 
23 5 

in (J fission probably reflects some characteristics of 

the emission mechanism of tritons which are different from 

239 

those of LEA. The Pu results have large errors due to 
poor statistics hence no definite conclusions can be drawn 
about the yield variation. 

One may therefore conclude that channel effects 
manifest themselves in the LEA yield and that the use of 
the LEA to probe the propertied of the fissioning nucleus at 
the saddle point have proved fruitful. 
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The method adoped in the present trajectory calcula- 
tions differs from most other calculations in that the 
initial parameters took a distribution of values i.e. no 
restriction was imposed on their initial values. The initial 
values (trajectories) which produced correct final values 
were selected by the method of assigning ’weights' P-sing 
experimental distributions and then selecting the relatively 
more 'weighted' trajectories. The initial distributions were 
obtained using these selected trajectories. 

The significant aspect of the results obtained is that 
they have shown that by and large it is possible to obtain the 
correct final distributions and most of the final correlations 
without imposing any restrictions on the variation of the initial 
parameters. Previous calculations have used either only low 
or only high initial energies. The authors had concluded that 
in order to produce the correct final distributions it was 
necessary to restrict the variation of the initial parameters. 

The initial energy distributions obtained in the present 
calculations are broad and show possible peaking in the low 
and high energy regions indicating that there is no apriori 
reason for choosing either only low or only high initial 
energies — both of these can produce the connect final 
distributions. 

Of the approximations used to simplify such calcula- 
tions, the point-charge approximation has time and again been 
held responsible for the inability in explaining or 
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reproducing some experimental observations (Boneh et al. , 

1967; Katase, 1968). It has been suggested that the inclusion 

of deformation effects would resolve these questions. An 

< 

attempt to include simple deformation effects in the present 
calculations has indicated that there is no significant 
change in the initial distributions. 

The experimentally observed amticorrelation between 
Ep and which is a unique characteristic of ternary fission, 
has been the crux of determining the validity of the choice of 
initial energies. Boneh et al. (1967) had concluded that in 
order to reproduce the correct anticorrelation it was necessary 
to use high initial energies whereas Choudhury (1976) has 
obtained the same result using low initial energies. The 
present results indicate that in order to obtain the correct 
anticorrelation it is necessary to restrict the variation in 
interfragment distance without regard to the initial energies 
being high or low. When the anticorrelation condition was 
used as a criterion in assiging 'weights', to the trajectories , 
the initial distributions remained unaltered indicating 
conclusively that the initial distributions are in fact broad. 

The idea that correlations between initial parameters 
do exist and that they effect the final correlations has often been 
been mooted (eg. Boneh et al. , 1967)* however, the few atempts 
that have been made in this direction, notably by Boneh et al. , 
have not resulted in any conclusive result on these correlations. 
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An attempt in this direction has been made in the present 
calculations and one result that has been obtained is that 
the initial alpha and fragment energies are positively 
correlated. The final negative correlation between these 
two quantities seemed to necessitate a correlation between 
the point of emission of the alpha particle and the inter - 
fragmont distance, however, it was not possible to obtain 
any conclusive result on this correlation. 
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